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1 Introduction 



The purpose of this thesis is the computation of the observables for heavy-light mesons, mesons 
with a light and a heavy valence quark, in QCD. There are various methods available to accomplish 
this task. Hence, the hrst question is which method is chosen for the computations. The degrees of 
freedom which should be used give the first constraint on the methods in question. In this attempt 
quarks and gluons should be used, the elementary degrees of freedom of QCD. This choice reduces 
the number of methods which can be used. In particular QCD Sum Rules (QSR) are chosen to do 
the job. 

Those few sentences raise several questions. There are two which are of special interests. How are 
the properties of resonances calculated in QCD the theory of strong interactions ? What are QCD 
Sum Rules ? The answers to these questions directly show how the calculations performed in this 
thesis work. Based on those calculations conclusions can be drawn. 

The classical computations based on the QCD Lagrangian have been perturbative ones. They 
where based on a perturbative solution of the QCD Lagrangian where the solution is expanded in 
the coupling constant a s and every order can be calculated via the Feynman diagram technique. 
Characteristic computations based on perturbation theory have been calculations of running cou- 
pling constants, running masses or of charges and deep inelastic scattering. Perturbation theory 
was also used to check whether QCD is the theory of strong interactions or not. 
Unfortunately the properties of mesons can not be calculated via perturbation theory. The per- 
turbative solution of the QCD Lagrangian is useful in the case where it can be truncated after a 
certain order. There are infinitely many orders of the solution in a s . A calculation of every order 
would be impossible. Therefore the coupling constant a s has to be small if the truncation should 
be valid. However, the coupling constant is momentum dependent. The coupling constant is small 
for momenta larger than Aqcd and large for momenta smaller than Aqcd- Thus, the perturbative 
solution is expected to be valid for large momenta and invalid for small momenta. This behavior 
is called asymptotic freedom because the interaction tends to zero if the energy goes to infinity. 
Asymptotic freedom is the property that allows for all perturbative calculations in QCD. 
In the time when the tests and measurements based on perturbation theory have been performed 
a new phenomenon was recognized. Free quarks or gluons have never been found. This was called 
the confinement of quarks. To our knowledge its impossible to produce free quarks or gluons. 
They are always confined in a bound state. Though it is widely believed that at high energies 
QCD matter melts and reveals quasi free quarks and gluons. This would be something like a liquid 
where the molecules are quarks and gluons. The confinement has up to today not been explained 
and is possibly not explainable by perturbative physics, but it is most probably close connected to 
hadron physics. 

The physics around us takes place in the domain of small momenta whereas the physics which 
takes place at high momenta could only be found in the early universe. There are artificial sources 
where the physics at high momentum scales can be found like accelerators. Hence, if the QCD 
physics at the energy scales at which we live should be described a non-perturbative solution of 
the QCD Lagrangian is needed. Unfortunately there is only one method which is supposed to be 
capable to solve the theory in every domain. It is lattice QCD, but this method needs a lot of 
computer power and is not very instructive about the physics which stands behind the calculations. 
Therefore, a less direct method is used, but a method with considerably more physical insight than 
lattice QCD. 

The dream of everybody who works in this field would be to find a formula as it is given for the 
hydrogen atom. Here every bound state has a main quantum number n which characterizes the 
mass of the bound state. Something like that is aimed for by modern researchers. The fine struc- 
ture of the problem is completely neglected because the raw structure is already that difficult. One 
method to represent the states are spectral functions where every bound state is identified with a 
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peak in a function of the energy. The properties of the bound states characterize the peaks. Thus, 
if the properties are known the peaks can be calculated or vice versa. The mathematical object 
which contains the spectral function of a certain bound states are correlators. This has been ex- 
plored in the K/"alcn-Lehmann spectral representation. In the domain with time-like momentum 
the correlator is given by the spectral function of the system. In a perturbative calculation bound 
states are not generated. 

However, there is a possibility to estimate the spectral function taking non-perturbative effects 
of the QCD Lagrangian into account. One exploits the fact that the correlator in the space like 
domain can be calculated by the so called operator product expansion (OPE). Via a dispersion 
relation time-like and space-like momenta can be connected with the domain where the momentum 
squared is negative. A phenomcnological ansatz for the spectral function is made and the proper- 
ties of this ansatz are fitted in order to reproduce the OPE. In a certain way the spectral function 
is then calculated by the scheme. QCD enters the scheme through the OPE which is the quantity 
where theoretical calculations enter the game. The so called Wilson coefficients of the OPE can 
be calculated using the QCD Lagrangian. Hence, this scheme has two parts a phenomenological 
part which is the spectral function and a theoretical part which is the OPE. This scheme is called 
a QCD Sum Rule. It has two features that provide the non-perturbative information. Through 
the OPE non-perturbative information enters the calculation and the dispersion relation connects 
the spectral function with the non-perturbative information. 
The thesis is structured in four parts: 

• In the first part, sections 2, 3, 4 and 5, the OPE and its calculation are introduced with a 
special glance on some subtle but important features of the OPE. 

• In the second part, section 6, dispersions relations are discussed. 

• The third part, sections 7, 8 and 9, introduce correlators, spectral functions, QCD Sum Rules 
and the Borel Sum Rules. 

• The fourth part, section 10, QCD Sum Rule calculations of heavy- light systems are presented. 
In this part novel calculations are discussed. In addition previous QCD Sum Rule calculations 
are reviewed. In these examples the Wilson coefficients are calculated in QCD and in the 
heavy quark effective theory (HQET). 

2 Regularization and renormalization 

In this section important terminologies and techniques in the field of regularization and renormal- 
ization are introduced. In the first part terminologies are discussed and in the second examples of 
some common techniques are presented. The results of the examples will be used in subsequent 
sections. 

2.1 The terminology of superficial divergence 

Let D be the superficial degree of divergence: an amplitude with the degree D diverges like: 




(1) 



Some examples are given here: 



D = 




(2) 
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This means: all amplitudes with D > are divergent, it is important to notice that also the 
amplitude with D=0 is divergent. Such amplitudes are logarithmically divergent. Examples for 
divergent diagrams are loop diagrams. A simple example occurs in the real scalar <f) 4 theory where 
the tadpole diagram is divergent: 




In this case D=2. The external legs have been amputated in the explicit expression. The number 
of loop integrals in each theory is infinite but fortunately there is a class of quantum field theo- 
ries in which the number of elementary diagrams which are divergent is limited. In this class the 
treatment of these divergences is simpler than in a class of theories with infinitely many divergent 
diagrams. An example for a diagram which is finite and contains a loop is 

qi q'i 

Here D=-4, the integral is not divergent. The external legs are again amputated and the external 
momenta are set to 0. 

2.2 Regularization and renormalization of the tadpole diagram in A0 4 - 
theory 

The Lagrangian of the theory is given by 

* = \ (Mf - A ^ ( 7 ) 

There is only one diagram to order A 1 (see equation (5)). This diagram is the simplest divergent 
diagram that exists. Divergent diagrams can be regularized. These regularized expressions pa- 
rameterize the divergence of the diagram. After the regularization it is possible to renormalize 
the diagram. Renormalization is a process in which the infinities are absorbed from observable 
quantities which have a relation to the amplitude, in order to eliminate them from the results. 
To achieve that a reference momentum is chosen at which the observable quantity is measured. 
This point is called the renormalization point. Taking the difference of the quantity given by the 
diagram at the renormalization point and another point results in an expression which is given by 
the value of the observable at the renormalization point and some momentum dependent function. 
At this stage the infinities cancel and the theory is defined at an arbitrary renormalization point. 
The full analytic expression for the diagram (5) is 

; ( _ a) f _dk^_i i_ = (~A) f dk* 1 



q 2 — m 2 J (2-7r) 4 k 2 — m 2 q 2 — m 2 (q 2 — m 2 ) 2 J (27r) 4 q 
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Ptotai = k- k + q- q = Q 



Figure 1: Momentum conservation at the vertex. 



The momentum conservation at the vertex in (5) shows that the external momenta do not con- 
tribute to the internal momentum (see figure 1). The diagram is factorized in the contribution of 
two free propagators and the loop. In the following only the integral is kept. This integral which 
represents the loop is divergent, its superficial degree of divergence is D=l. There exist various 
methods to regularize and renormalize divergent diagrams. Every method has advantages and dis- 
advantages therefore the three most common are applied to the tadpole diagram for a comparison 
and an overview. 

An explicit integration of the loop integral is needed. One method to accomplish the integration 
is the Wick rotation. Here the integration in the four dimensional Minkowski space is transformed 
to an integration in a four dimensional Euclidean space. The only thing employed is the residue 
theorem. 

The integration over the contour in figure 2 is zero. In the limit R — oo only the integration 




Figure 2: Contour of integration in the zero component of (8). The red dots are mass poles shifted 
by e away from the real axis. If the limit R — ► oo is taken the integration over the bows vanish. 



along the real and the imaginary axis contribute. Therefore the integrations must be equal except 
of their signs. This enables the replacement of the integral along the real, I\, axis with the integral 
along the imaginary axes, I^. The integral I2 is not affected from the poles sketched in figure (2) 
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and furthermore the replacement transforms the integral from Minkowski to Euclidean space. This 
transformation is further investigated in section 4. The whole equation can also be merged in a 
graphic picture. The integration contour of I\ can be rotated into the integration contour of I2 
without crossing the mass poles. Hence the name Wick rotation is explained as a simple rule for 
the procedure just outlined. 

After the Wick rotation the integral is given by the following expression: 

The pole in the integrand has vanished and the integration runs in Euclidean space, this is sym- 
bolized by the index E. The transformation law is simply 1° = il E , I = l E , the minus sign in front 
of the integration was explained above. 

Thus the integral is computed and the comparison of the renormalization methods can be started. 
1. The cutoff method: 



This method is the most simple one, application of the Feynman rules results in an integration 
from to 00. This integral is divergent and has therefore to be regularized. Regularization 
means a parameterization of the divergence or in other words a method to make the integral 
finite with a rule that tells how the original integral can be derived from the regularized one. 
Thus one way to achieve that is by replacing the upper limit by A. The result is: 



dk 4 



1 



(2tt) 4 k 2 



16tt 2 



(A 2 -m 2 ln 


"A 2 


hi 


) 




TO 2 







forA^>m 



16tt 2 
A 2 -m 2 ln 



m 2 ln[k E 
A 2 ' 



in 



T 



-if (A). (10) 



The original divergent result is retained in the limit A — > 00. This is a parameterization 
of the tadpoles divergence. Through this procedure a parameterization of the divergence is 
obtained, the integral is regularized. Without the cutoff the integral is infinite. 
Renormalization is a method to get rid of the unphysical parameter A, which is infinite and 
therefore unphysical. The question is which effect this divergence has. As we know from 
measurements the propagation of a particle from a to b is not infinitely likely. This gives us 
the strong hint that the divergence should be absorbed into something inside the amplitude. 
If we neglect all other contributions to the propagator except of the tadpole we can form a 
geometrical series, that looks like: 



D(p 2 )=i J d 4 x e <p " I "<O|r{^(O)0(a:)}|O) 



p2 _ m 2 



l-i/(A) 



p2 _ m 2 



-*/(A) 



p2 _ m 2 



+ ... = 



p A — m z 1 _ 
i 



-*/(A) 
p 2 — m 2 

i 



p 2 -m 2 +if(A) p 2 - m 2 R 

This shows that the mass can absorb the divergence which we parameterized by the cutoff. 
The combination m 2 — if (A) is called the renormalized mass rriR. The massparameter that 
we called m is the naked mass and has to be divergent in order to absorb the divergent 
—if (A). If we assume the naked mass to be divergent, this subtraction of divergences elimi- 
nates the problem of infinities. Thus renormalization of the diagram is achieved. 
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The renormalization just performed has to be treated with care. The loop integral is in- 
dependent of the external momentum p^. This means the renormalization of the mass by 
means of the tadpole is a shift of the mass it does not produce a momentum dependent 
mass, the renormalize mass stays constant. The first momentum dependent loop emerges in 
second order A. Therefore no physical measurable effect exists that proves whether the mass 
is shifted by the tadpole or not. 

The consequence of this analysis is a simplification of our calculations! All loop diagrams 
in which the loops do not depend on the external momentum p can be dropped from the 
perturbation expansion of the two point function. They produce no momentum dependence 
of the mass but the mass is a parameter which enters the theory out of measurements. This 
means it is fixed if all contributions beyond first order A can be neglected. 

2. The Pauli-Villars method: 



Another method to regularize the tadpole diagram is the introduction of regulators. Where 
the regulators arc just copies of the original particel's propagator with a very big "mass". 
The starting point is the loop integral 

d\A 1 r Jhi 1 r ^1.4 



(2tt) 4 k 2 
with the conditions 

and 



i n f dk* _J f dk_ 

-m 2 1 J (2 7 r) 4 fc 2 -M 1 2 2 J (2tt 



1 



) 4 fc 2 - Mf 



D 1 + D 2 = 1 



D X M\ +D 2 M% = m 2 . 



(12) 

(13) 
(14) 



The additional propagators are the regulators with "heavy masses" M,. After the integration 
these conditions arc directly needed. 



16tt 2 



A 



; ln 



A 2 1 
— + 1 



DiA 2 + D X M\ In 



A 2 



M 2 



D 2 A 2 + D 2 Mf ln 



A 2 



M| 



16tt 2 



A 2 - L>iA 2 - D 2 A 2 - to 2 ln 



"A 2 " 




"A 2 " 




" A 2 " 


m 2 


+ L>iM 1 2 ln 


Ml. 


+ D 2 M% ln 





Together with (13) A 2 - L>iA 2 - D 2 A 2 = and 



DiM 2 ln 



A2_ 

Ml 



D 2 Ml ln 



Ml 



(DiMf + D 2 M%) ln 



"A 2 " 




"A 2 " 


= m 2 ln 


_M 2 _ 


_M 2 _ 



the result for the regularized loop integral is 

—i 



16tt 2 



m 2 ln 



M 2 



(15) 



(16) 



(17) 



In comparison with (10) the power term A 2 does not appear. Therefore this term is regular- 
ization dependend. 

However, the renormalization procedure shown in the paragraph about the cut-off method 
holds here too. The physical quantity, the mass does not change, because the shift due to 
the loop is unobservable. 

The phenomenon just observed is of general nature and can be comprised as follows. Different 
rcgularization procedures can result in different expressions for the same diagram, but the 
renormalized observable quantities are not allowed to change if the regularization procedure 
is changed. Unless this condition is satisfied the theory used would make no sense at all. 
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3. The method of dimensional regularization: 



Another method to regularize the divergence in (8) is to extend the integral in 4 dimen- 
sions to an integral in n dimensions. The extended integral is divergent for even n with 
n > 4. The loop integral in n dimensions is given by 



dk- 4 _„ i 1 ,4-»r(i-f) 



2 r 
16^2 



(2tt)™^ k 2 -m 2 (4tt)9^ ( m 2)i-f 

2 /m 2 \ 
7 - 1 - ln(47r) + In — + 0(4 - n) 

- — n V M / 



(18) 



where 7 is the Euler-Mascheroni constant and /i 4 ~™ is introduced to keep the dimension of 
the expression constant with respect to n. The integration can again be performed with the 
Wick rotation. Thus, the third method to regularize the integral gave another expression 
that differs from the results of method one and method two. All methods up to now gave 
expressions that where different from each other. Hence, the expressions are regularization 
dependend. 

The dimensional regularized loop diagram has to be renormalized. Therefore it has to be 
subtracted. Two schemes to do that are very common the MS and the MS scheme. In the 
MS scheme just the pole term —jzr^ is subtracted while in the MS scheme — 7 is 

subtracted. Therefore the renormalized integrals are scheme dependent. This example shows 
that not only regularization but also renormalization can be scheme dependent. 
Again the physical quantities are not affected from such scheme dependences. The tadpole 
diagram in dimensional regularization (18) has to be summed up as shown in the cut-off 
method before it is subtracted. The scheme dependence does not show up in the renormalized 
mass. It is just shifted from the bare to the physical value. 

The methods shown here are convenient for the calculation of one loop diagrams. A criteria which 
constrains the choice of the renormalization technique is if the symmetries of the theory in which 
the calculations are performed are still valid after renormalization. 



3 The operator product expansion (OPE) 

This section provides a basic introduction of the operator product expansion. The expansion is 
derived in two ways where the first is the physical introduction and the second one the mathematical 
introduction. After that issues which are important in the application of the OPE are considered. 
Finally questions how the OPE should be truncated and how it can be simplified are addressed. 

3.1 Features and their Derivation 

The purpose of this section is to establish a first overview of the OPE and topics connected to it. 
The most important features of this expansion are explained and derived in the following sequence. 

1. The Operator Identity 

In 1969 Wilson proposed a Non-Lagrangian model of current algebra [7]. The basis of his 
model is the assumption that the singular part as x —* y of the product A(x)B(y) of two 
operators is given by a sum over local operators 

A(x)B(y)^J2 C n B ( x -y)°n(y) (19) 

n 
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where C^ B (x — y) are singular c- number functions and O n are operators. This relation is 
derived by using the formalism of the path-integral. This access to the OPE gives valuable 
physical insight about the OPE, in particular on the separation of long and short distance 
fluctuations. The remarkable thing about the OPE is that it is an operator relation; this 
means, in applying it to any matrix element (/3\A(x)B(x)\a) , the same functions C^ B (x — y) 
are received for all states |a) and \(3). 

2. Derivation in the Operator Formalism 

The OPE can also be derived in the operator formalism of quantum field theory. This 
derivation deals with the singularities of the operator product A(x)B{y) in the limit x — ► y. 
These singularities have to be separated from each other. This operation produces a series 
in which the summands are sorted by the strength of their singularities which coincides with 
their importance. The series derived in this scheme is identical with the OPE derived by 
using the path-integral formalism. 

3. Important and Unimportant Summands 

An expansion with infinitely many summands that can not be resummed is useless if it 
can not be truncated. Using scale arguments it can be shown that the OPE can be trun- 
cated. 

Dimensional analysis suggests that C^ B (x — y) behaves for x — > y like the power do — dA — ds 
of x — y, where di is the dimension of the operator i in powers of mass or momentum. Since di 
increases as we add more fields or derivatives to an operator O, the strength of the singularity 
of C^ B (x — y) decreases for operators C of increasing complexity. 

The decrease of the singularity in equation (19) with operators 0(y) of increasing complexity 
is the reason that justifies the truncation of the OPE. Generally speaking the OPE is useful 
in drawing conclusions about the behavior of the product A(x)B(y) for x — > y. 

4. Translation to Perturbation Theory 

For application purposes the OPE is translated into the language of Feynman diagrams. 
The Wilson coefficients are displayed as a series of graphs that contain common symbols and 
symbols representing the connection to the matrix elements. While the matrix elements are 
simply multiplied with the expression for the corresponding Wilson coefficient without any 
instruction how they should be calculated. Special attention is payed to the fact that high 
and low frequencies have to be separated. This leads to the introduction of a renormalization 
point fi. 

3.2 The Operator Identity 

In this section a generalized version of the Wilson expansion is derived. The derivation is close to 
[25]. For this purpose, let us consider a Greens function for local operators A\(xi), A 2 (x 2 ), etc. 
whose arguments approach a point x, as well as other local operators B\(yi), B 2 (y 2 ), etc. with 
fixed arguments: 

(0\T{A 1 (x 1 ),A 2 (x 2 ), B 2 (y 2 ), ...}|0) 

= J \U l , z d(t)i{z)}a 1 (x 1 )a 2 (x 2 )...b 1 (y 1 )b 2 {y 2 )...ex P (iI[((>\) (20) 

where the lower-case letters a and b indicate replacement of the field operators, the A's and B's, 
with c-number fields cj>. The operators A\(xi), A 2 (x 2 ), etc. form the operator product which 
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should be expanded, while the operators B\ (j/i ) , B 2 (y 2 ) , etc. represent external states. The c- 
numbcr fields <f) are the field operators of the considered theory. 

Now surround the point x with a sphere B(R) of radius R which is much larger than the separations 
among the x,x\,x 2 , etc. but much smaller than the separations between x, yi,y 2> etc. (see figure 
3). Since the action is local, it may be written as 



/ = f d A z^(z) + I d A z^{z). 

JzeB(R) JzgB(R) 



(21) 



Equation (20) can be put in the form 

(0\T{A 1 (x 1 ),A 2 (x 2 ), ...B 1 (y 1 ),B 2 (y 2 ), ...}|0) = 

/ \H-i,z$B(R)d<t>i(z)\ bi(yi)b 2 (y 2 )...exp I i / d 4 z£f(z) 

J \ JziB(R) 



/ [ni, z eB(R)d(f>i{z)] a 1 (x 1 )a 2 (x 2 )...exp I i / d i z^f(z) 

J \ JzeB(R) 



(22) 



in which the only influence that the fields inside and outside the sphere have to each other is 
given by the boundary condition at the surface of the sphere. The two fields have to be connected 

e 2 



X^XXXX X'X^XX XX 


KX XX XXX X XX X 
.r KaXaX 







ei 



Figure 3: Partitioning of K 2 into a Ball of radius R around the point x and the remainder of the 
space. This picture has to be treated with care because in the 2 dimensional Minkowskispace a 
sphere is not a sphere. This point is going to be elaborated in section (4). 



continuously to each other. Thus, the boundary conditions are that the field or the n-th derivative 
of the field is continuous on the surface of the sphere. The fields inside and outside the sphere have 
the same value for a certain point u on the sphere's surface. Except of this, the dynamics of the 
fields are completely independent from each other. They could even be described by two different 
theories. Hence the amplitude has been split into two domains. 

This fact is going to be exploited to simplify the expansion that was just performed. The integral 
over the fields inside the sphere are fully determined by the values and derivatives of the fields on the 
surface of the sphere, which in turn may be expressed in terms of the fields and their derivatives 
extrapolated from outside the sphere to the interior point x. Hence, after this conversion the 
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integral over the sphere will not depend on specific coordinates of the spheres surface but just on 
its radius R. 

In the simple example that is shown in figure (4) the extrapolation is quite easy. The field 




V V V V V v V v v v V 



ei 



Figure 4: A simple example of a field that can be split up as shown in figure (3). The field is 
constant outside the sphere but it has a gradient inside the sphere. 



is a constant outside the sphere so the extrapolation to the point x will be the same constant. 
This constant has already been the boundary condition on the spheres surface. In conclusion the 
boundary conditions stay the same. 

The integral over the fields inside the sphere contains the part of the amplitude where short 
distances are involved while the integral over the fields outside the sphere contains the part of the 
amplitude where long distances are involved. In addition to this separation the two domains do 
not share any variables. 

If we express this integral as a series in products of o(x) and the c- number fields <f> and their 
derivatives extrapolated to x, the coefficients can only be functions C^ 1 '" 42 '"' (x\ — x, x 2 — x, ...) of 
the coordinate differences and R. Since the points yi, y 2 , etc. are all far outside the sphere, the 
sphere has no effect in the limit R -C \x — so in this limit equation (22) becomes 

(0\T{A 1 (x 1 ),A 2 (x 2 ), ....Bi(yi), B 2 (y 2 ), ...}|0) - 

J b 1 {y 1 )b 2 {y 2 )...exp(i J d 4 zJ?(z) 

x ^2 C n ,: - r : -x,x 2 - x, ...)o n (x) 

n 

= E C n 1,A2 '-(*i -x,x 2 - x, ...)(T{O n (x), B^yt), B 2 (y 2 ), ...}) (23) 



for x\, x 2 , etc. all approaching x, where 0(x) is the quantum-mechanical Heisenberg-picture 
operator corresponding to o n (x). The operators that occur in this expression depend solely on x. 
Therefore, they are called local operators, a non-local operator would depend on x and further 
coordinates. 

In comparison with (22) the integral over the fields inside the sphere has changed drastically. 
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It has separated into summands which are in turn separated in two factors. The first one is a 
coefficient function while the second one contains c-numbcr fields and their derivatives. While the 
integral over the fields outside the sphere stayed nearly the same. The operator O n that stems 
from the integral over the fields inside the sphere connects the long and short distance part of the 
amplitude. 

This yields by Fourier transforming (23) with respect to the y variables and multiplying with 
appropriate coefficient functions 

(0\T{A 1 (x 1 ),A 2 (x 2 ), ...}\a) -> ]T Ct' M '-{xi -x,x 2 - x, ...){/3\O n (x)\a) (24) 

n 

for arbitrary states \a) and (/3|, of course the states correspond to the B operators. Because this 
applies for arbitrary states, it is the operator product expansion in a generalized version: 

T{Ai(x\), A 2 {x 2 ), ...} -» Cn uA2 '-(xi -x,x 2 - x, ...)O n (x). (25) 

n 

It is very important to understand that the OPE arises only in the limit R <C \x — yi\. This is a 
special and not the general case. That means the applications of the OPE are limited to matrix 
elements where you have the distinct separation between the distances \x — Xi\ and \x — yi\ of the 
form 

\x - Xi\ < \x - yi\ . (26) 



3.3 Derivation in the operator formalism 

The derivation in section 3.2 provides a very intuitive introduction of the OPE. Many features 
of the expansion can simply be read off the proof. Some of them where already claimed. Still 
missing is the scheme that orders the sum appearing in the OPE. Also the analytic properties 
of the coefficient functions have not been analyzed. This is done in the present section via an 
alternative derivation of the OPE. This section is close to [9] . 
A product of two operators is considered 

P(x ll ,^)=T{A 1 (x lt +^),A 2 (x lt -^)}. (27) 

The vector £ M is split into spherical coordinates: 

p = VW, ^ = y ( 28 ) 

Henceforth P is considered as a function of x^, r/^ and p : 

P{x^^) = P{x^,p). (29) 
If P diverges for p — » we define an operator 

1 {x^) = Y™ P{x ^ p) (30) 

in the weak limit, the limit performed for the matrix elements and not for the norm of the operator. 
This limit should be finite. Hence, C\(p) has to compensate the divergence of P(x fl7 r]^, p) at p = 0. 
Therefore, C\{p) has to have a singularity at p = by itself. C\{p) is a suitable function. The 
singularity of C\ is restricted by the condition that the result be finite and different from zero. 
Under the assumption that there are some matrix elements of P 

{^\P{x' V 'p)W) (31) 
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which near p — are as singular as or more singular than all other matrix elements of P a suitable 
Ci can be found. The mathematical term for this requirement is 

($|P(a: M , w )|qr) 

exists for any tp 1 x^ and r]^. Choosing any of these most singular matrix elements, equation (31), 
as the function C\ the operator (30) should have finite and non-vanishing matrix elements. 
By (27) and (30) the operator 0\ is local in the sense that 

[A j (x),O 1 (y,r ] )} ± =0 
[O 1 (x, V ),O 1 (y, V ')} ± =0 (33) 

for {x — y) 2 < with commutators or anticommutators taken appropriately. 

In many cases 0\ turns out to be a multiple of the identity. In perturbation theory it seems to be 
a general rule that the most singular part of a matrix element is given by 1 

($|P(a;,rj,p)|*) « (0\P(x, v , p)|0)<$|tt) (34) 

provided the vacuum expectation value does not vanish 

(0\P(x, V ,p)\0)^0. (35) 

This leads to the determination of 0\ 

lim <*l p (*.^)l*> = <^> (36) 

or 

Oi = 1. (37) 

If the vacuum expectation value vanishes 

(P(x, V ,p))=0 (38) 

(30) may lead to a suitable definition, but we will see that in general there are more local operators 
which should be associated with P. 

We now turn to the second problem of analyzing the singularities of the operator product. To this 
end the remainder is introduced 

r 1 (x, V ,p) = P{ *^P) _ 0l[Xjr]) (39) 

which vanishes in the limit 

lim ri(x,77, p) = (40) 

because of (30). In order to get some information about the singularities of P near £ = we 
multiply (39) by C\(p) and solve for P 

P(x,r ] ,p) = C 1 (p)0 1 (x,r ] ) + P2(x,ri,p) (41) 



1 There may be several equally singular matrix elements. This possibility is not considered here. 
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with 

P 2 (x,r 1 ,p) = C 1 (p)r 1 (x,7],p). (42) 

Here C\ — > oo and r\ — > for p — > 0. Hence nothing can be said in general about the limit of P 2 . 
If 

lim P 2 = (43) 

p^O 

(42) already gives complete information on the singularities of P near £ = 

P(z, = Ci (p) Or (x, r/) + P 2 (x, v , p) (44) 

with liirip^o P2 {x,7],p) — for 77 constant. Hence the term C\0\ carries the total singularity of 
P. 

If P 2 diverges the product P has additional singularities near £ = 0. In that case its possible to 
repeat the procedure for P 2 that was just applied to P with 

P 2 (x, V ,p)=P(x, n, P ) - d (p)Oi (a;, r,) ■ (45) 

This time the most singular matrix element C 2 of P 2 is chosen to form 

2 (x, V ) = hm — . (46) 

Hence, more information on the singularities of P is acquired. Introducing the remainder 

P 2 (x,T],p) 

r 2 (x,t],p)= c ^ 2 (x,r7) (47) 



with 



lim = r 2 (x,r],p) = 0. (48) 

p^O 



An Operator analogous to P 2 can be introduced 

P 3 = r 2 C 2 . (49) 

In conclusion it is obtained 

P 2 = C 2 2 + P 3 . (50) 

Inserting (50) into (41) it is obtained 

P(x, n, P ) = Ci(p)Oi (x, 77) + C 2 (p)0 2 (z, 77) + P 3 (x, 77, p). (51) 
C 2 is less singular than C\ 

lim^g =0. (52) 
p^o 61 (p) 

This procedure has refined our analysis of the singularities of P, (52) follows from 

lim C2[P) lim iteiM _ 1im W 2 (x, Vl pM = 

P^o(0|P 2 (x,77,p)|V)/C 2 (p) P^o(0|P 2 (x,77,p)|V)/C 2 (p) 1 ! 
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because lim p ^o »"2(x, ??, p) = and lim p ^o P cQp) P ^ = C>2(x, 77) 7^ for at least one matrix element. 
This process can be iterated until a P n is derived which vanishes at £ = 0. It is an assumption that 
the process terminates after a finite number of steps. The result of the iteration is the expansion 

P(x,r],p) = Ci{p)Oi(x,ri) + ... + C n (p)O n (x,r])+R{x,r],p) (54) 

with the following properties 

hm = (55) 

lim C n (p) = co or lim C n {p) = c 7^ (56) 

p^O p^O 

lim R(x, r), p) = 0. (57) 

p^O 



The operators 3 are given by 



n. = i im P T,iJiC a O a 
J p^o Cj 



and satisfy the locality conditions (33). 

The choice of the matrix elements that are used as the Cj is not unique, in turn the operators Oi 
are not unique. This arbitrariness can be reduced by requiring that the Operators Oj are linearly 
independent. The transformation that has to be performed to achieve this can always be arranged 
without changing the general features of the functions Cj . An expansion with linearly independent 
Operators Oi and coefficient functions that have the features listed in (55)-(57) is called a standard 
expansion. 

This alternative derivation of the OPE gives insight in the analytic properties of the coefficient 
functions. Obviously the analysis of singularities in the operator product P(x, rj, p) at p — provide 
an alternative approach to the OPE, in addition to the separation of the amplitude in short and 
long distance fluctuations. 



3.4 Important and unimportant summands 

The operator products P that are investigated have a certain dimension d. Wilson's assumption 
about the OPE leads to the statement that only operators O n with dimension < d are relevant in 
the OPE. This statement is derived by scale arguments. The chain of arguments is shown in the 
following [9]. The starting point is the assumption that a scale invariant theory is investigated. 
This means the theory is invariant under the transformation: 

x — ► sx. (59) 

In the quantum mechanical sense this means that there is a family of unitary transformations U (s) 
with the property: 

U^is^^Uis) = s d ^(j>(sx) (60) 

for any local field operator of the theory. The real number d(<p) is called the dimension of the 
operator (f>. Using this transformation law Wilson determined the singularities of the coefficients 
in the following way. Applying the transformation U (s) to the expansion 

00 

A,(x + 0A 2 (x - = E C i^)OM (61) 
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we get 



U- 1 (s)A 1 (x + Z)A 2 (x-$U(s) = J2u-\s)C j (Z)O j (x)U(s) 

3 = 1 

CO 

U- 1 (s)A 1 (x + 0A 2 (x -£)U(s) = J2 Cj {^U-^Oj (x)U(s) 



3 = 1 

OO 

„d{A 1 )+d{A 2 ) 



Ai{ax + s£)A 2 {sx -s£)=Y^ Cj(Z)s d(0i) Oj(sx). (62) 

3 = 1 

Expanding AiA 2 on the left hand side yields 

oo 

s ^)+^)J2C j (sOO j (sx) = J2C j (Os d ^ ) O j (sx). (63) 

3 3 = 1 

Since the Oj are linearly independent the following equations are valid 

- s d(A 1 )+d(A 2 )-d(O j ) c .^ ^ 

as the scaling law of the coefficients. The exponent of s is called the dimension of Cj and we have 
the result that in each term of the expansion the dimension of Cj and the dimension of Oj must 
add up to the total dimension of the left hand side 

d(Cj) + d(Oj) = d(A t ) + d(A 2 ) = d{A l A 2 ). (65) 

The dimension of Cj indicates the behavior for £ — > 0: 

CM) = s- d{C,) Cm (66) 

or 

C j {x) = \x\- d{C * ) Cj{^). (67) 



As a consequence the Cj(x) are singular or ^ in the limit x — > only if d(Oj) < d(A\A 2 ). Hence 
the singular and non-vanishing terms are provided by operators Oj of dimension less or equal to 
the dimension of the product A\A 2 . 

All the terms belonging to operators with d(Oj) > d(A\A 2 ) vanish in the limit x — > 0. The 

conclusion is that all Operators with d(Oj) > d{A\A 2 ) are not relevant for the OPE. 

The simple power-counting argument above is modified by renormalization effects; the expansion 

(19) must be formulated in terms of operators renormalized at some scale n, and then fj, appears 

along with x — y in the coefficient function C^ B {x — y) . An important example are asymptotically 

free field theories, where C^ B (x — y) does behave like the power do — dj± — d B of x — y suggested 

by dimensional analysis only up to a power of ln(x — y) 2 . 

The corresponding statement in momentum space is that for k — > oo, 



d 4 xe- ik " x »A(x)B(0) -» J2 V n B {k)O n {0) (68) 

and correspondingly 

' d 4 xe- ik " x »T{A(x)B(0)} -> ^ C* B (k)O n (0) (69) 
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where V^ B (k) and C^ B (k) are functions of k^ that for large k 2 decrease more and more rapidly 
for more and more complicated terms in the series. 

A short remark concerning the validity of the arguments given above has to be made. The basis of 
all arguments have been the equations (59) and (60) . Unfortunately dimensional arguments do not 
work in every theory, in some theories they are not valid. In all realistic theories scale invariance is 
broken as an example any mass term breaks the scale invariance. Furthermore it can happen that 
the dimensionality d can not be defined in a usefull way. Hence in such theories scale arguments 
can not be used in order to determine if the OPE can be truncated. 

However, dimensional arguments work very well in perturbation theory. For the exact solution 
the situation may be quite different. This can be demonstrated in the Thirring model. There 
Wilson's hypothesis holds in perturbation theory with the conventional dimensions d{A\A2) = 
d{A\) + d(A 2 ), but this is no longer true for the exact solution of the Thirring model [8]. Hence, 
the conclusion is that the arguments above work in perturbation theory and perturbation theory 
will be used in the remainder of this thesis. Therefore, the arguments can be used. 



3.5 Translation to perturbation theory 

The preceding sections about the OPE have been very general. Calculations in quantum field theory 
are mainly done via perturbative techniques. Therefore this section introduces the perturbative 
technique to calculate the OPE of a Greens function. 

The first step is to define the operators in the OPE. The operators are given by combinations of 
the field operators appearing in the Lagrangian of the theory. In standard expansions (see section 
3.3) the operators are linearly independent. That means the operators in standard expansions are 
given by all possible combinations of field operators of the theory under consideration, that are 
linearly independent. In this and the following sections only standard expansions are used. These 
arguments fix the operators. 

The second step is the calculation of the coefficients in Wilson's expansion, the so called Wilson 
coefficients. The recipe was already given in section 3.3. There the Wilson coefficients have 
been identified with the coefficients of matrix elements in the operator product expansion. The 
singularities of these coefficient functions determined the operator they belonged to. This method 
has two weaknesses. The Wilson coefficients have not been determined uniquely. There can be 
several different Wilson coefficients which have singularities of the same strength. In perturbation 
theory the strength of the singularities of these Wilson coefficients is often not determinable. An 
alternative method is urgently needed. 

The clue that provides this method is that the OPE is an operator relation and the operators O n 
are already known because we use a standard expansion. If an OPE is sandwiched between two 
external states 

(a\P\0) = (a\^2C n O n \0) =^2(a\C n O n \0) =^C n (a\O n \0) (70) 

n n n 

the Wilson coefficient C„ is unaffected by the external states. This is the feature that is exploited 
in the calculation of the C„. There are states that filter out certain coefficients when the OPE is 
sandwiched between them. The matrix elements of the operators O n with these states are all zero 
except of one or a few of them. In this section the perfect case in which just one matrix element 
is not zero is stressed. Suppose the matrix element with Oj is not 

( a \p\p) = cmoM ^ Cj = 

{a\Oj\P} 

This is a perturbatively calculable expression but there is a restriction that requires a modification 
of the calculations usually made in perturbation theory. In the first derivation of the OPE, section 



19 



3.2, it was shown that the coefficients C n stem from a path-integral that is restricted to a Ball 
with radius R in four dimensional space time. Therefore only short distance fluctuations play a 
role in the calculation of the coefficients. This is the point that allows the Wilson coefficients 
in many theories to be calculated via perturbative techniques. The dynamics which takes place 
at short distances happens at large momenta. In asymptotically free theories this is the domain 
where perturbative calculations are valid. The impact of this observation is that no small momenta 
contributions occur in the Wilson coefficients. The scale \x that distinguishes high and low momenta 
is defined by additional conditions. 
The concrete rules that follow are 

1. The external momenta must be bigger than /i, q? > [i. 

2. The virtual momenta in the loops of involved diagrams must be bigger than \x. 

Otherwise the OPE is not applicable to the problem. The method just introduced is called the plane 
wave method. With these rules all Wilson coefficients are calculable, although the calculations are 
sometimes very difficult. There exist alternative methods but the plane wave method is very 
valuable because it provides a lot of physical insight in comparison with other methods. 

3.6 OPE in the real scalar 4 theory 

In the last sections the groundwork concerning the OPE was done. This work should be substan- 
tiated in quantitative calculations. Additionally there appear complications like the question of 
the application in theories with broken symmetry These points are treated in this section. The 
calculations should be as easy as possible. For this reason the most simple theory is used, the real 
scalar <f) 4 theory. The Lagrangian density is given by 

* = \ W \<<t> 2 (72) 

where the dimension of <j> is energy. This theory has an intrinsic symmetry, it is a simple reflection 
symmetry 

JSf (0) = J2f (-0) . (73) 
The Feynman rules are given in table 1. As an example the operator product (fxp is expanded in a 

vertex : — i\o 
propagator : 2 _ \ 



Table 1: Feynmanrulcs in the <j> -theory in the phase where the symmetry is not broken. The 
vertex is a four vertex. 

standard expansion. The set of linear independent Operators is given by 

. (74) 

Operators with an odd number of 0's can be excluded because they do not respect the reflection 
symmetry. The important terms in the OPE are the terms with d < 2 (see section 3.4). With 
these informations the OPE is determined to be 

T {0(g), <£(<?)} = Cil + C 2 «X0)<X0). (75) 
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This is the OPE for the propagator in the interacting cf) 4 theory. The OPE is going to be derived 
in first order in the coupling constant A . It is convenient to start with the left hand side of the 
OPE (75). The vacuum expectation value is given by 

D(q) = (0\T{^q)^(q)} |0) = -i J d 4 xe^ x » (O|T{0(x)0(O)}|O) 

Q— 



1 + 



m. 



(76) 



o 



(see section 2). This Greens function is the sum of the two terms in the OPE (75). The OPE can 
be obtained with the rules given in section 3.5. First, the requirement that q is large is exploited. 
The propagators are expanded in 
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In the propagator we keep terms, up to order in order to be consistent with the second term in 
(76) 
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Secondly, the loop integration has to be separated in a high and a low frequency part. This 
corresponds to a separation in short and long distances. Hence a scale [i is introduced which 
separates high and low frequencies, high frequencies > \i > low frequencies 
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(79) 



This is the OPE of the propagator, although the components, Wilson coefficients and operators 
remain to be identified. The expansion could have been done in an alternative way by using the 
plane wave method, discussed in section 3.5, to calculate the Wilson coefficients C\ and Ci. The 
plane wave method exploits the operator nature of the OPE. If the OPE is sandwiched between 
appropriate states a certain summand in the expansion is singled out. Based on that the Wilson 
coefficient of the summand can be calculated, but the calculations can also be comprised into 
diagrammatic techniques. If this is done the Wilson coefficients are expressed as diagrams. The 
analytic result corresponding to the diagrams are the Wilson coefficients. The states for C\ would 
be (0| and |0) and for C 2 (p\ and \p). The calculation of C 2 is simple, the corresponding diagram 
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is given by 



(80) 

where the cross denotes the contact with the operator. The calculation is simple, its just the 
product of two propagators and the four vertex. This product is taken in the limit where the 
external momentum is large. Hence, a Taylor expansion in is performed and the lowest order 
terms are kept. After all a symmetry factor of 2 has to be taken into account. The result turns 
out to be C 2 = ^r- 

The matrix element (O|</>(O)0(O)|O) can be calculated via known techniques. It is just the Greens 
function (q\(p(0)^>(0)\q) in leading order Ao where the external propagators and the four vertex are 
amputated. The amputated part is just the Wilson coefficient Ci. An example for the contractions 
is given by 

^ J dM g\mmHpW(p)Hp )Hm- (si) 

The corresponding diagram is given by 



(82) 

where the Wilson coefficient C2 has to be amputated. The cross symbol denotes a special vertex, 

it is the operator <f)(0)(f)(0). The result is just — , the expression for the operator in lowest 
order. An alternative way to derive this result would be to calculate the tadpole diagram in <f> 4 
theory with the extension that the frequencies large than /i have to be cut off. Hence the third 
summand in (79) is 6*2(0 |0</>|O), with respect to (75) the remainder has to be C\\. Thus the 
following identifications can now be placed: 
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C 2 = P, (85) 
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A remarkable point is that a running mass arises in the OPE in first order in Ao- The mass is 
running with respect to the scale [i which separates high and low frequencies. In the literature 
H is referred to as the renormalization point (see [10]). Without this separation a running mass 
does not appear to first order in Ao, it would be constant without the separation in high and 
low frequencies. The ultraviolet cutoff M and the naked mass mo arc unobservable quantities 
and should not appear in the definition of the running mass. Therefore, the running mass is 
renormalized in order to eliminate the unobservable quantities. 

The Wilson coefficient C\ was calculated to first order in Ao- The first order result can be used to 
calculate a part of the higher order corrections to C\. In higher orders, diagramms appear which 
consist of n first order results which are ordered one after another 

+ Q + — Q— Q— + .... (87) 



Hence, a series analogous to the geometrical series is derived 
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is identified as the mass of the particle. Obviously, the divergent parameters m and M still appear 
in the expression for the mass. These parameters have to be eliminated in order to make the theory 
well defined. A new scale /x is introduced 



m 2 (/x) = m 2 , 
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Hence 



m 2 {n) = m 2 (/Lt ) 



32tt 2 



(95) 



Thus all divergent parameters have been eliminated. The mass is renormalized and everything is 
well defined. Hence, (83-86) is the OPE of the propagator in the 4 theory to leading order in the 
coupling constant A. 

This simple example tells a lot about the OPE. The matrix elements are not purely non-perturbative. 
They receive contributions from perturbation theory. The OPE reproduces the perturbative result 
to a certain order in the coupling constant A and in 

3.6.1 OPE in the phase with broken symmetry 

A second look at the Lagrangian (72) of the 4 theory exhibits an interesting feature of the theory 
if the mass is exchanged by a mass parameter with arbitrary sign 



2 = \ (dtf 
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The Hamilton operator is given by 
H = J d 3 x 
The potential is given by 

y(0)=±^V + ^*. (98) 
Depending on the sign of the ?7 2 -term in (96), the potential has the shape of figure 5. In the case 



1 





Figure 5: The potential of the Hamilton operator in 4 theory, depending on the sign of r] 2 . On 
the left hand side the case of a positive sign is shown while on the right hand side the case of a 
negative sign is shown. 



of a negative sign in front of r\ the symmetry of the system is spontaneously broken. The order 
parameter which determines this is (O|0|O). If (O|0|O) ^ the symmetry is broken. In the case of 
broken symmetry it is reasonable to expand the theory around the minima of the potential 

= 0o + o-(aj). 0o = ±v = ±J^-ir}. (99) 
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unshiftcd theory 


vertex : — i\ 
propagator : p2 _' m2 


shifted theory 


vertices: —i\ ,—i^J^-ri 
propagator: p ,_' 2T)2 



Table 2: Comparison of the Feynmanrules of (/> 4 -theory in the phase where the symmetry is not 
broken with the rules in the phase where the symmetry is broken. An additional vertex arises in 
the phase with broken symmetry. 



Where </>o is given by the position of the minima in the <fr space. The Lagrangian S£ after this 
transformation is given by 

* = \ (^) 2 - \f^W 3 ^V- (100) 

The Feynman rules in this version of the theory differ from the ones in the theory where the 
symmetry is not broken. A comparison is given in table 2. Where — i\q is a four vertex and 



— iyj ^rj is a three vertex and is of order VAq- If the OPE instead of the usual techniques are used 
to calculate the amplitudes in the phase with broken symmetry the use of the new Feynmanrules 
can be avoided. There are two possibilities to incorporate the nature of the states in the calculation 
of the amplitudes. The first one is to shift the theory to the minima in the <f> space which results 
in new Feynmanrules. The second is to use the OPE and incorporate the nature of the state not 
in the Feynmanrules but in the matrix elements. This is especially helpful if for some reason the 
shift which leads to new Feynmanrules is not possible. In the following it is proven, that the OPE 
reproduces the amplitudes in the phase of the theory with broken symmetry. 
The OPE can directly be taken from section 3.6. Obviously, only m 2 , has to be replaced by — rj 2 
and the matrix elements change. Again, it is emphasized that the theory is not shifted. Hence, the 
Wilson coefficients do not change in comparison with section 3.6 except for the mass parameter. 
Thus the Wilson coefficients stay the same, only the matrix elements change. 
In the following a first naive test is performed if the OPE reproduces the amplitudes in the phase 
where the symmetry is broken. This test will fail, but in the next section the reason for this failure 
will be clarified. 

The method to calculate the matrix elements in the OPE in cj) 4 theory is to use the classical value 
for </> which is given by 4>q. In this case 

(0^(0)0(0)10)=^^) =^ 2 . (101) 

This is the value for (0|(/>(0)(/>(0)|0) at the scale \x at which "m 2 hysical (ji) = —2i] 2 , as it is already 
stated in equation (100) in summand 2. 

The treatise just performed leads to the following OPE in the phase with broken symmetry 



i J dxe~ ix ^ (O|0(a;)0(O)|O) =1 + ^-^1 
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This expression should reproduce perturbation theory but how is the propagator in the phase with 
broken symmetry calculated ? In order to keep the calculations as simple as possible the shifted 
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theory should be used. The diagrams that contribute to first order in Aq are 




(103) 



Only the last diagram contributes to the propagator the other ones vanish after renormalization 
(see section 2). In the limit of big momenta q the expression for the propagator is 



and it does not coincide with (101) even after setting [i = q and ^ = m v h ys . This result is of 
course shocking. The question appears if everything that has been stated about the OPE in the 
phase with broken symmetry is wrong. Fortunately it is not wrong but there is something missing 
in the derivation of the OPE. In the next section the riddle about the missing component will be 
solved. 

3.6.2 Comparison of OPE with perturbation theory in the phase with broken sym- 



The question that has to be answered is if the OPE reproduces the amplitudes in the phase 
with broken symmetry to a given order in Aq. For this purpose the propagator is expanded via 
perturbation theory in the phase with broken symmetry and via the OPE. The results should 
coincide to a given order in Ao but as shown in the preceding section they seem not to do it. The 
missing parts of the OPE are derived in the following. 

As shown in section 3.6 the normalization point is crucial in the definition of the OPE. The missing 
agreement between the OPE and perturbation theory in the phase with broken symmetry could 
also stem from an incomplete implementation of the renormalization group flux. In the version of 
the theory with reflection symmetry the only scale that is important is the mass in the Lagrangian 
(72) which naturally defines the scale (i which is inherent to the theory but it is not necessary 
to use the physical mass of the particle as the renormalization point. An other point is that the 
physical mass in the theory with broken symmetry is not given by the mass parameter 77 in the 
Lagrangian (96) but by — 2?7 which is the natural scale in the phase with broken symmetry. 
These complications did not show up in the theory with unbroken symmetry where there exists no 
difference in the mass scale that the Wilson coefficients are calculated with and the mass scale that 
is set by the physical mass. Thus the flux should be less important in the theory with reflection 
symmetry than in the theory where this symmetry is broken. 

In simple words the Wilson coefficients are calculated in a version of the theory that is defined on 
an other scale than the version of the theory in which the propagator (104) is calculated. Thus 
the point at which the Wilson coefficients are defined have to be changed to the point where the 
propagator (104) is defined. Unless this is done the effective theory (100) and the OPE will not 
coincide. 

In order to make the OPE independent from the normalization point \i the flux of the Wilson 
coefficients and all parameters has to be constructed. The method used here to achieve this is the 




(104) 



metry 
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method of renormalization group improvement. The OPE is improved by summing up the leading 
logarithms using the recipe that renormalization group theory proposes. 

The first step is the renormalization of the coupling constant Ao . Those computations are standard 
and can be taken from books like [22]. The renormalization group equations are given by 



„ d _ \{p) 



(105) 



H—Ci = o 2 



%TT 2 



(106) 



(107) 



d W X o A2 (A*) 



(108) 



All of these equations except of (105) are derivatives of the first order expressions. The derivation 
of (105) leads to new insight concerning the renormalization group flux of the OPE summands. 
The derivation can be found in [11] and is not shown here. The solutions of the equations (105-108) 
are 



AG") 



ln 



(mo) 



(109) 



m 2 (n) = 



a(m) 



m 2 (/z ) 



(110) 



C2 (M) = 



A(mo) 



C 2 (no) 



(111) 



Ci0i) = Ci(/io)-2 



m 2 (^o)C l 2 (^o) 
A(Mo) 



A(mo) 



(112) 



This is the improved OPE in the theory with reflection symmetry corresponding to the theory 
with a positive value of the squared mass parameter in the Lagrangian (96). It has the remarkable 
ability that it reproduces both versions of the propagator in the <j) A theory. If fi = mp it reproduces 
the propagator in the theory with reflection symmetry and if /j = — 2toq the propagator in the 
theory where this symmetry is broken. This assertion is proved in the following. 
The starting conditions are chosen to be 



^ / 1 m 2 (q) 

Q 1 



(113) 



C2(no = q) 



2<7 4 



(114) 
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then the coefficients C\, Ci contain no logarithms. The explicit expression for the OPE with those 
conditions is 



D(q) = (0|C 1 (m) + C 2 (m)^ 2 |0) = ^ 



1 



A(g) 
A( M ) 



2m 2 (/i) 



A( M ) 



A( M )<</> 2 ( M ))-m 2 (/z) 



(115) 



With these results the validity of the OPE in </> 4 theory can be tested. 

1. Coincidence of the improved OPE with the propagator in the theory with reflection symme- 
try: p? = m 2 

At this scale the matrix element given in (86) vanishes which leaves the OPE as 



[ A(«) " 


1 

3 


[ A(g) 1 


_A(m)_ 




A(m)_ 



1 I 

which coincides with the propagator in the theory with higher symmetry to first order in A. 
2. Coincidence of the improved OPE with the propagator in the theory with broken symmetry: 



4 -m 2 (m) + 0(A 2 ), 



(116) 



p? = —2m 2 = m 2 



phys 



At this scale the matrix element is given by (101) which leaves the OPE as 



D{q) 



q 2 q 4 



-2m 2 (m ph y S )) ( 2 
1 1 



q2 q4 m phys 



\ A(g) ' 


3 


[ A(g) 1 


.A(m p hj, s )_ 




X(m p hys) _ 


[ A(g) " 


2 
3 


\ A(g) 1 


_ \(m p hy S ) _ 




_A(TOph ys )_ 



1 



m 



phys 



3A , 

1 H o l n 

16tt 2 



Tflpfoyt 



+ G(A 2 ), 
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which reveals the coincidence with the theory in the phase with broken symmetry. 



^ 2 = to 2 



/i 2 — TOp/^s 



Figure 6: An illustration concerning the distribution of normalization points. 



The conclusion of this section is that the change of matrix elements, claimed in section 3.6.1, is 
not sufficient to change between the phases of the theory. In addition the normalization point has 
also to be changed. 

The ability of the OPE to describe processes in the phase with broken symmetry simply by changing 
the matrix elements and the normalization point in comparison with the OPE in the phase where 
the symmetry is not broken is the reason why it is employed in modern applications. 
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In the early days of the OPE the question arose whether the OPE is an approximation to the 
theory in the phase with broken symmetry or if it is exact in reproducing the amplitudes [10]. The 
source of this problem was an approximation made by Shifman, Vainshtein and Zakharov in their 
original publications to QCD Sum Rules [1]. This approximation is discussed in the next section. 

3.7 OPE in QCD - Wilsonian OPE versus Practical OPE 

The OPE has been explained in general and in terms of an example in the real scalar </> 4 theory but 
every theory has its own peculiarities concerning the OPE. In the remainder of this thesis QCD is 
the theory in which calculations are performed. 

Concerning the OPE in QCD the OPE exhibits big advantages. QCD also has a broken symmetry, 
the chiral symmetry. Unfortunately, the theory is much more difficult than the 4 theory. In 
the phase with broken symmetry an additional phenomenon arises, the confinement of quarks. 
This means, that free quarks do not exist in that phase. Only bound states of quarks and gluons 
have been measured, the so called hadrons. No one knows how the transition between the phase 
with unbroken and broken symmetry has to be described. Thus a shift as it has been done in 
section 3.6.1 does not seem to be possible. Therefore, no Feynmanrules for amplitudes on the basis 
of quarks and gluons in the phase with broken symmetry are known. There exist effective field 
theories, which describe the physics in that phase on the basis of hadrons as degrees of freedom. 
Moreover, there exist other methods like lattice QCD which seeks for a numerical solution of the 
QCD Lagrangian. An additionally method to describe the physics of QCD in the broken phase is 
given by the OPE. 

The OPE is for example an alternative to chiral perturbation theory. As it was stated before the 
Wilson coefficients can be calculated in the phase with and without broken symmetry. Hence, 
one can use the QCD Lagrangian, where the quark masses are current quark masses, to compute 
the Wilson coefficients. A repetition of the calculations performed in the <j> 4 theory would lead to 
amplitudes valid in the phase with broken symmetry. In this expansion quarks and gluons are the 
basic degrees of freedom. This is the way the OPE is used in actual applications. 
There are two approximations that simplify the calculation of the OPE in QCD considerably 

1. The perturbative parts of the matrix elements can be absorbed in the Wilson coefficients 
of the unit operator 1. Thus, in this approximation the matrix elements are purely non- 
perturbative (see [2]). 

2. The renormalization group flow of the Wilson coefficients can be neglected. 

In the propagator example of 4> A theory the perturbative part of the 4> 2 {0) matrix element has been 

_ 2 2 

In • Thus in the first approximation, this term is included in C\ and must be subtracted 
from the matrix element. 

Those two approximations are applied to all OPE computations in QCD. They date back to the 

founders of the QCD Sum Rule method [1] and have lead to missunderstandings [10]. 

The exact OPE is denoted Wilsonian OPE and the approximative one is called practical OPE or 

SVZ expansion/OPE. In the remainder of this thesis the practical OPE is used. 

One generally can describe the OPE as a method to approximate the physics in the nonperturbative 

regime, the phase with broken symmetry, with the field theory that is valid in the perturbative 

regime. The big advantages in this case are that the OPE deals with current quarks and not with 

constituent quarks. 

Every advantage has its price. In the OPE case it is the introduction of the matrix elements which 
in QCD are called condensates. Up to today there exist no method which is based on first principles 
to calculate these matrix elements. Chiral perturbation theory and QCD Sum Rules can be used to 
get an approximation for the condensates, but they are based on data from expirements and not on 
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the QCD Lagrangian. Hence, the condensates are determined through measurements. Nevertheless 
the OPE is heavily used in nonperturbative physics. In the next section the calculations of the 
Wilson coefficients are shown for several examples. 



4 Calculation of Wilson coefficients in QCD 

In the last section the OPE was introduced and it was showed that the Wilson coefficients are 
the part of the OPE that has to be calculated if QCD amplitudes are concerned. The plane-wave 
method was introduced as the technique for the computations, but the details of the calculations 
have not been explained. Therefore, this section is devoted to the calculations of Wilson coefficients 
via the plane-wave method. This method is the most simple one can think of, it directly employs 
the operator nature of the OPE. Moreover, the first calculations of Wilson coefficients where done 
with this method [I]. As it was quoted the practical OPE is used, thus the renormalization group 
flow of the coefficients is neglected. 



4.1 Equal mass case : the flavor conserving scalar current correlator 

Here and in further investigations products of currents are concerned. The currents are used as the 
interpolating field for a hadron. Unfortunately, this choice is not unique. As the approximative 
current for a hadron all currents have to be considered which have the same number as the hadron. 
Such ambiguities can have different implications. For example the ambiguity between a four quark 
and a two quark current has a physical implication. The quarks in the current are the valence 
quark. Hence, such an ambiguity is directly linked to questions concerning the quark model. Non- 
physical ambiguities are for example given by the number of derivatives which occur in the current. 
It is possible to write down n quark currents with and without derivatives which have the same 
quantum numbers. In this thesis mesons are analyzed. Thus, the currents have to have a even 
number of quark fields, but here only two quark currents are addressed. Moreover, currents with 
the simplest structure are used. This means for example currents with the smallest number of 
derivatives. The most simple current of this type is the scalar current: 

j(x) = q(x)q(x) j f (x) = q(x)q(x). (118) 

Hence the OPE is: 

oo 

j(x)jH0) = Y,CnO n . (119) 

n=l 

The current product has energy dimension 6, hence section (3.4) states that only operators with 
dimension < 6 appear in the OPE. Those operators are given by [1] 

1 dimension = 
mqq dimension = 4 
G C ^ U G C ^ dimension = 4 
A c gG^„ dimension = 6 
qTiqqT2q dimension = 6 
f^G^G^G^ dimensioned. (120) 

To be precise there are more operators with dimension < 6 than those just quoted, but they can 
be reduced to those quoted above. An example of such operators are operators with a derivative. 
Derivatives raise the dimension of an operator. Hence, a derivative of the mqq still appears in the 
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OPE but it can be reduced to the mqq Operator. The calculations of Wilson coefficients gets more 
complex if more complex operator products are concerned, for a given operator product it gets 
more complex the more complex the operator is to which the Wilson coefficient belongs. Therefore, 
a simple example for the calculation of Wilson coefficients is the Wilson coefficient that belongs to 
the operator mqq, the quark condensate. According to the plane-wave method this coefficient can 
be filtered out by sandwiching the OPE between quark states: 



(p\j(x)j*(0)\p) = ( P \ ]T C n O n \ P ) = £>|C n O„|p> - Yl C n(p\On\ P ) = 
n—1 n—1 n—1 

C« q ( P \O qq \p) = C qq (p\qq\p) = C qq u{p)u{p). (121) 

The current correlator is transformed to momentum space 

i J #xe* x {Q\T{j(x)?(QW) = * J d 4 xe^(0\T(q(x)q(x)q(0)q(0))\0). (122) 
After that the plane-wave method is applied for the calculation of the qq Wilsoncoefficient 

i J d 4 xe^(p\T(j(x)jHO))\p) = i J d 4 xe^(p\T(q(x)q(x)q(0)q(0))\p). (123) 
The coefficient is derived to lowest order in as- Therefore Wick's theorem is applied: 

l J d 4 xe^(Sx)q(xMO)qmP)+^ J d 4 xe^(p\q( x)q(x)q(0)q( 0)\p) 




i(p + 4 + to) _ i(f — d + to) 

= S \ p+q) 2_ m 2 U ^ S ) + ™(P, S ) (p _ g) 2_ m 2 M fa S )' ( 124 ) 

In the diagrams dashed lines represent currents and full lines represent fermions. Hence, an ex- 
pression for the left hand side of (121) is 

i> + d + to ii — d + m 

The goal is to reduce those expressions to a form where one can identify the quark term uu. The 
external particles are free particles p 2 = to 2 . Hence, the Dirac equation u{p){f — m) = or 
{f — m)u(p) = can be applied 

TO+| + TO TO-jf + TO 

~ u{p \p + q y-m 2 u{p) - u{p \p- q y-m 2 u{p) = 

2to + d _ 2m — d 
~ u (p)i \o — — o u (p) - u ip)l ^5 o u {p)- (I 26 ) 

(p + q) - m {p - q) - m 2 
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What about the ^-term ? The treatment of this term is a bit more difficult ! 



Tfl I I 

u(p)jfu(p) = u(p) — tfu(p) = 7^u(p) {mtf + 0m} u(p) = ^u(p) {ft + #»} u{p) 



2m 



1 



—P^qMp) {'f'f + 7 V} «(p) = —p^q v u{p)2g' lu u{p) = 
Inserting the last expression in (126) results after a trivial conversion in 
1 2m 2 + p^q„ 1 2m 2 



u(p)u(p) 



[mu{p, s)u(p, s)} 



p^q^ 



m 2 (p + q) 2 — m 2 

The Wilsoncocfficicnt is determined to be: 

1 / 2m 2 + 
to 2 \ (p + q) 2 
2m 2 +pV qfl 



m 2 (p — q) 2 — m 2 



mu(p, s)u(p, s)} . 



(127) 



(128) 



p^q^ 



2m 2 



p^q^ 



m z 



p 2 + q 2 + 2p )1 q^ 



+ 



(p - q) 2 

2m 2 



p^q^ 



p 2 + q 2 - 2p^f 



m" 



(129) 



From the derivation of the OPE we know that the Wilsoncoefficients can only depend on q 2 . The 
dependence of this Wilsoncoefficient is now analyzed. The coefficient takes the form: 



m* 



2m 2 + p»q^ 2m 2 - p"g„ 



(p + q) 2 - 
1 / 2to 2 



p^q» 



(p - q) 2 - 
2m 2 



p^q^ 



m 2 \ q 2 + 2p f _ t qt 1 



2P^ 



(130) 



The only p dependence enters the coefficient through the scalar products. This dependence should 
not occur in the coefficient after all that we know. The dependence stems from the fact that we 
only take the effect of a single plane- wave into account. This wave has a fixed four momentum 
but there is no physical circumstance that marks out a special momentum vector. This means we 
have to take all possible vectors into account. This is achieved by transforming the scalar product 
to a space with a Euclidean metric. In this space the scalarproduct is averaged over the Euclidean 
angle leaving the coefficient without any p dependence. 



4.1.1 Transformation to 4-dimensional space with Euclidean metric 

The most convenient way to perform the upcoming integration is to transform the problem from 
the Minkowski-Space to a space with a Euclidean metric. This is done by decomposing the -1 in 
the metric tensor of the Minkowski-Space 

(131) 

9i±v — n n _1 n n n _1 n (132) 



The — i 2 is subscribed to the component of the 4 vectors in which the tensor is sandwiched. The 
effect of this procedure is that the component which was the component becomes imaginary and 
will from now on be called the 4 component. The advantage of the effort is that in the space in 
which we have transformed the problem the metric is Euclidean. This enables is to perform the 
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integration mentioned above. 
The transformation laws are: 



9n» 

P 2 - 



- -1 4 

2 2 
Pe = -P 



Pa — ► Pea = ipo- (133) 
The action of the transformation can be illustrated in two dimensions. The metric transforms as: 

(134) 





' 1 




" 1 


" 


9»u = 


-1 







1 



while the vectors transform as: 



Po 
Pi 



PE„ 



Pi 

ipo 



(135) 



(136) 



Following the common conventions there is a zero component in the Minkowski, but not in the 
Euclidean space. The zero component in the Minkowski space becomes the last component in the 
Euclidean space. 

An interesting effect concerning this transformation is the transformation of the sphere, p 2 — const.. 
The sphere in the Minkowski space is a hyperbola while the sphere in a Euclidean space is a ball. 
The transformation (133) transforms these two objects into each other (see figure 7). In four 




Pi 




Figure 7: The sphere in the 2 dimensional Minkowskispace and the sphere after applying the 
transformation 133. 



dimensions the transformation causes exactly the same. Thus, from this example its plausible that 
the transformation enables the use of common techniques for the averaging process. 

4.1.2 Averaging over the 4-dimensional Euclidean angle 

In the preceding section, the method used to average over all possible orientations of p^ and 
was explained. In this section an example and the Wilson coefficient are computed. 
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n 


<(jv«T> 
(IpIMP 


2 


1 


4 


1 


6 


fi4 



Table 3: Averaging the scalar product. 



The four dimensional Euclidean measure is: 

dp% = \pe\ 3 sin 2 (e 2 )sin{e 1 )de 2 de 1 d(j)dp E . (137) 

The solid angle in 4 dimensions is: 

dn 4 = sin 2 (e 2 )sin(9 1 )de 2 de 1 d(t). (138) 
While the surface of the unit sphere is given by 

/>77 />7T f-27T 

/ / dn A = 2tt 2 . (139) 

J 8-2=0 J6i=0 J<t>=0 

If the average over the scalar product of two vectors pe^Me.^ in the 4 dimensional Euclidean space 
should be calculated it is very convenient to define the 4 axis in the direction of the vector qE lfJ ,- 
Then the scalar product reads: 

p E ^qE.^=cos(e 2 )\p E \\qE\- (140) 

Therefore averaging over arbitrary orientations of these two vectors means to perform the following 
integration: 

Se 2 =o jgUo J^l P E, ' l <lE,ndn 4 
fe 2 =o /<?i=o f<t>=o d ^ 

\PE\ \q E \ ' ' 



2tt 2 



S 2 =0 J6»i=0 



cos(e 2 )sin z (92)sin{e 1 )de 2 de 1 d(l) = 0. (141) 



This trivial example leads to the formula for averaging {pE^q E ) n '- 

(\ I I \\ n r n r 71 /*27r 
< (p^) n >= U ' / / / co* n (02)«n 2 (0 2 )«n(0i)d02d0id0. (142) 

2tt J6» 2 =o Je 1= o J4>=o 

For odd n the result vanishes. For even n some results are shown in table 3. This example 
directly shows how (130) can be averaged over all possible orientations of p^ and q^. After the 
transformation to Euclidean space the procedure just outlined is repeated. The only difference is 
the integrand in (141) it is replaced by (130). This integration is not trivial but possible, algebraic 
computer tools can perform it easily. The result is 

Hi-,)( 1 + 3 .) 

rrr 1 + v 



where v — y 1 — ^f-. Some authors have an additional power of m in the coefficient [21] because 
they take qq and not mqq as the condensate. 

An important remark still has to be stated. The integration that has to be performed during 
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the averaging process could here be done in closed form. In general this is not possible. Thus in 
such cases the integrand has to be expanded in the scalar product. The term in the expansion 
can all be averaged as it was shown in (141). Therefore an expansion of the Wilson coefficient 
is calculable. At this point the problems begin. The summands are in principle given by powers 
of if to is small compared to q 2 the expansion can be truncated. The first summand is often 
a satisfying approximation. In the case of big masses to this can not be done, each summand is 
important. Hence, the series has to be resummed. The resummation is a difficult process. During 
the investigations made for this thesis no systematic method was found to achieve the resummation. 
However, the authors of [1] somehow solved this problem. 

The renormalization group flow of the coefficient is neglected as it was explained in section 3.7. 
Thus the calculations are finished and can be comprised as the sequence: 

1. Apply the plane- wave method with appropriate external states. 

2. Bring all free spinors to the right of the expression and identify them as placeholders for the 
condensates. 

3. Eliminate the external momenta, either by setting their square equal to the corresponding 
masses or by using an appropriate average method. 

4.2 Equal mass case : the flavor conserving pseudoscalar current corre- 
lator 

In this case everything is analogous to the scalar current correlator, except of a trivial point 
concerning 75 matrices. The current is given by 

j(x) = q{x) (i<y 5 ) q(x) j\x) = q(x) (175) q{x) (144) 

where the i is necessary in order to make the current self adjoint. Again the coefficient is calculated 
to lowest order oP s . The application of Wick's theorem results in 

i J d 4 xe tq ^Sx)(i l5 )q(^0)(i l5 )qm P )+i J d 4 xe lqx {p\g{ x) (t 75 )g(x)g(o) (i 75 )g( 0) \p) 




ft + 4 + m _ ft — d + m 

= -u(p) (l 75 ) - ^ 2 _ m2 (175) U(p) - U(p) («75) ( p _ g )2 _ m 2 (*T>0 ( 145 ) 

Step 1 is this time a bit longer because the 75 matrices have to be accounted for. Their action is 
the change of the sign in front the slashed quantities. The expression which determines the Wilson 
coefficient is altered into 

— t — 4 + TO — i) + d + TO 

aM (p' .)?-„■ ■"" + "ft" ft,-, )?-„■ ■"'"• (146) 
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Now, using the same manipulations as for the scalar current the spinors can be transferred to the 
right 



Pul I , Pu<l l 



u{p)u{p) + — -?p —u{p)u{p) = 



(p + q) 2 — m 2 {p — q) 2 — m 

— o o P ^ mu(p)u(p) H k P ^ mu(p)u(p). (147) 

m 2 q 2 + 2p^ y ' y ' m 2 q 2 - 2p^ K ' K ' v ; 

In comparison with (128) the numerator of (147) has obviously changed. Hence, the expression to 
be averaged is given by 



1 -v»(f , 1 p^q 



m 2 q 2 + Ip^qV m 2 q 2 — 2p f ^q tJ - 



(148) 



The relevant integral is of the same type as discussed in section 4.1 and can be obtained in closed 
form 

c-^^W- 1 - ( 149 ) 

4.3 Unequal mass case : the flavor changing scalar current correlator 

From the OPE point of view mesonic systems consisting of two quarks with unequal masses are 
more complicated. In the case of scalar mesons such systems are interpolated by the current 

j(x) = q 1 (x)q 2 (x) ]\x) = q 2 {x)q 1 {x). (150) 

In order to calculate a Wilson coefficient of such systems the sequence just quoted has to be 
passed trough as it was done in the case of equal mass systems. The principle of the calculation 
remains unchanged, but the steps are more complicated. The vacuum polarization transformed to 
momentum space is: 

i J d*xe^{Q\T{j(x)?{G))\Q) = 1 J d i xe lc ' x (0\T(q 1 (x)q 2 (x)q 2 (0)q 1 (0))\0). (151) 

For the calculation of the mqq Wilson coefficient the plane-wave method is applied 

i J d 4 xe^( Pl \TU(x)jH0))\Pi) = * J d 4 xe i ^( Pl \T(q 1 (x)q 2 (x)q 2 (0)q 1 (0))\p 1 ). (152) 

The plane- wave method discussed in section 4.1 must be altered a bit to account for the two 
different types of quark fields. Thus, two quark condensates enter the calculation, one for each 
field. Correspondingly, the external states of the amplitude are either of type 1 or type 2. This is 
symbolized by the 1 at the external states. Again the Wilson coefficient is derived to lowest order 



; J d 4 xe^ (p5^ (x)q 2 (0) gl (0) \p, 1) 



q+p 

q 4 V 



i> + <t + m 2 

=-^( < + ; )2 _ m ^ )- («3) 
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Note that the exchange term does not contribute to this amplitude. Using the Dirac Equation 
u(p) (fl — mi ) = or u\ (p, s) (p — m\ ) =0 leads to 

mi + d + TO2 

-Ul(P) , , x 2 2 Mp)- ( 154 ) 

The j^-term is treated as in section 4.1. The amplitude is given by 

( • \2 2 — Mp)Mp)- (155) 

m l (P + 9) ~ TO 2 

In subsequent applications the limit mi — > is of interest. Therefore it will be calculated here 

1 ml + mim 2 + p^ _ . . . . . . 

o 9 9 — ;r u i\P) u i\P)- (156) 

mi g 2 - m 2 . + m 2 + 2,^ w w v ; 

The expression for the Wilson coefficient is expanded in a series 

1 mj+mim 2 +M M _ m 2 + m 1 m 2 + p^ 



mi (g 2 - ml) ! + "i+2«> mi (g 2 - m 2 ) 



1 -^^r + -- (157) 



In the limit mi — ► only the leading term survives, remember p^q^ = miy/q^qi T cos(6) because 
PE.iP l E = PuP 1 * P er definition 

_ m 2 + \ q \cos(8) _ 
q z - m 2 

After averaging over the four dimensional Euclidean angle 9 the term is given by 

(159) 

q z — m 2 

Hence the Wilson coefficient in this approximation is given by 

-TO 2 m 2 

C * = -qT^ml = mFV' (16 ° } 

4.4 Unequal mass case : the flavor changing pseudoscalar current cor- 
relator 

The calculation goes along the same lines as in sections 4.2 and 4.3. The computations are illus- 
trated here: 

j(x) = qi(x) (i-y 5 ) q 2 (x) j\x) = q 2 (x) (175) qi (x). (161) 
Hence, the Wilson coefficient is determined by 



i / d 4 xe^(p| 1 g 1 (c C )(i7 5 )g 2 (x)g 2 (0)(i75)gi(0)b) ] 



1 1 



q + p 

q 4 \ q 
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In comparison with (155) only the sign of the TOiTO 2 -term in the numerator is changed. Thus, the 
result in the approximation where mi = is given by (160) with m 2 replaced my — m 2 

4.5 Equal mass case : the flavor conserving vector current correlator 

Vector mesons where the valence quark and antiquark have the same flavor are extrapolated by 
the current: 

3n(x) = q( x hM x ) J't( x ) = q( x hM x )- ( 164 ) 

The standard plane-wave analysis discussed above leads to the expression: 

(p\Ux)jUo)\p) = 



d A xe^(p\q{x) liiq (x)q{0)luqmP) +i I <fxe«>* {p\q(x)^q(x)q(0) lu q(0)\p) 



, { i ■ m _ i> — d + m 

= "(P» {p + q)2 _ m 2 ^<v) + u{p) lv (p _ g)2 _ m2 7 M u(p). (165) 



The occurrence of open Dirac indices is new. Such indices did not occur in the calculations to scalar 
particles. Fortunately the expressions for vector particels can be factorized in Lorentz tensors and 
Lorentz scalar parts as: 

IV (i 2 ) = {9,, qj f-) n T ( g 2 ) - M^n L (g 2) . (i 66) 

Due to current conservation only the transversal part contributes to the amplitudes, the longitu- 
dinal part can be neglected. Therefore it is removed from the amplitudes. After this subtraction 
only the transversal part times a Lorentz tensor is left. The transversal part is a Lorentz scalar. In 
order to keep the calculations simple only the transversal part is calculated. In comparison with 
the scalar and pseudoscalar coefficients the additional work is given solely by the procedure which 
extracts the transversal part of the coefficient from the full expression of the coefficient. Three 
steps are necessary in order to accomplish the extraction. 

1. Contract the amplitude with g^, and average this expression over the four dimensional 
Euclidean angle. The result is 3Cg ?i T (q 2 ) — Cqq.L (<? 2 )- 



2. Contract the amplitude with ^r-, and average this expression over the four dimensional 



9m 9" 

Euclidean angle. The result is —Cq q .L (q 2 )- 
3. Subtract those two expressions and divide by 3. The result is Cq q: T (q 2 ) = Cq q (q 2 )- 

Some comments on the calculation of Wilson coefficients can also be found in [21]. As an example 
the Wilson coefficient for the quark condensate in the vector current correlator case is calculated 
to lowest order in the quark mass. 

1. Contraction of (165) with g^ v 

_, , {f + j)+m 2 (p-4)+4m 

u ( p \ P + q y- m i u{p) + u{p) (p- g )»-m* u{p) = 



u(p) — 
m 



pq m 2 +pq 



(p + q) 2 ~ m2 (p ~ <z) 2 



u(p). (167) 
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The expressions have to be averaged over the four dimensional Euclidean angle. This can be 
done directly with the expression in (167) or with each term of its Taylor expansion. The 
expansion is done in the scalar product pq, the result to lowest order in the mass is 



4m 



2. Contraction of (165) with 



u{p) f (p + + m) jf | u(p) $(p-tf + m 
-u(p) 



(p + q) 2 



q 2 (p - q) 2 



t(p) {2qp<j l + q 2 {-f + <j t + m) 



(p + q) 2 



u{p) 



l(p) (2qp^ + q 2 (-p-^ + m)) 



{p-qf 



u(p) 



m 



2(pq) 2 



(p + qf 



q 2 pq 2{pq) 



i 



q 2 pq 



(p-qf 



u{p) = 
u(p) = 
u{p). 



(168) 



Again the Taylor expansion is considered. In lowest order of the mass the longitudinal part 
of the Wilson coefficient is ^r 1 . 

3. The last step is simple and the result is 



C- m ^ 



(169) 



The solution quoted in (169) can be used for light quark systems. In the case of heavy quark 
systems all orders in he Taylor expansion have to be considered. Hence, the expansion has to be 
summed up or the average has to be taken direct without doing a Taylor expansion. The result 
where all orders in mass are included is: 

n (a 2 )- 2 ( 1 ~^( 2 + ^) (170) 

The examples just shown illustrates the plane-wave method for the computation of Wilson coeffi- 
cients. They have been arranged from simple to difficult but it does not end at the vectorparticle 
level. Tensor particels of arbitrary rank, can also be considered. With increasing rank the Dirac 
structure and the corresponding projector are more and more complex. Additional complications 
arise in loop diagrams and in diagrams involving gluons. Furthermore, there are ambiguities in 
the interpolating currents (see section 4.7). 

The plane-wave method is the most rudimentary method for the computations of Wilson coeffi- 
cients. Every method posses advantages and disadvantages. The advantages of the plane-wave 
method are that it can be used for the calculation of every Wilson coefficient without having to 
learn a great amount of new formalisms. The disadvantage is the complexity in the calculations 
which have to be performed. There exist methods where the calculations are significantly less com- 
plex. The most important one besides the plane-wave method is the fixed point gauge technique, 
also called background field method. There are many examples concerning the application of this 
method in [15]. 
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4.6 Diagrammatic representation of an OPE 



An OPE of the two point correlator in a meson channel can be graphically represented. This is 
illustrated below for the first two terms, in the OPE 



(p\j,(x)ji(o)\ P )=-..- 




■ ■ ■ {mqq). (171) 



This includes all terms in the OPE to lowest order in as- The term in front of the unit operator is 
referred to as the perturbative term. This terminology stems from the fact that the practical OPE 
(see section 3.7) is used where the Wilson coefficient of the unit operator is given by perturbative 
QCD. 

The term in front of the quark condensate is similar to the scattering diagram used in the calculation 
of the corresponding Wilson coefficient above. In fact the diagrammatic representation can be 
understood as stemming from the calculations done in the plane-wave method. The crosses which 
do not appear in the scattering diagram symbolize the contact of the Wilsoncoefficient with the 
quark condensate. 

Note that the connection between the scattering diagrams and those representing the Wilson 
coefficients is non-trivial due to the angle average in Euclidean space discussed above. 

4.7 Quark mass effects 

The last sections have been very general and mathematical but for applications of the OPE also 
physical boundary conditions have to be taken into account. A very important boundary condition 
is given by the mass of the quarks that are chosen to build up the interpolating currents. They 
are important for Wilson coefficients involving the gluon condensate. 

Suppose the Wilson coefficient for the gluon condensate is calculated for the case of a current 
which is build up of two quarks with different flavor. Then two masses mi and mi enter the 
Wilson coefficient. Depending on the meson which should be approximated by the current various 
limits for the masses are interesting. In the case of a heavy-light system the limit 777,2 — > is 
interesting. Where mi should be the mass of the heavy and m 2 the mass of the light quark. 
Surprisingly this limit exhibits a serious problem. The Wilson coefficient in this limit diverges! 
According to physical and mathematical arguments this divergence should not occur. Hence, 
further considerations are necessary. Methods to handle them where already starting to develop 
in the original publications to the OPE [1], became fully developed in [5] and where completed by 
[21]. 

Here a short introduction is given following [5]. Sandwiching the OPE between one gluon states 
gives a suprising result: 

Y,C n (k\O n \k) = C mi (k\m 1 q 1 q 1 \k) + C m2 (k\m 2 q 2 q 2 \k}+C G (k\G 2 \k} + .... (172) 

n 

Calculating the matrix element (k\mqq\k), that means writing down its OPE leads to 

(k\mqq\k) = ±-^(k\G 2 \k) + ... . (173) 

We see from this expression that the quark matrix element is of the same order in a s as Cq , since 
in lowest order the coefficient C m is of zeroth order in a s . Inserting (173) into (172) the following 
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expression is obtained 



J2c n (k\o n \k) 



U G T 777 °mi + 77: ^m 2 

12 7T iZ 7T 



(fc|G 2 |fc) + ... . 



(174) 



Therefore the coefficient of the gluon condensate differs from what it was expected to be. The 
physical version is given by: 



1 a s 

Cg, physical = Cg + "777 (Cmi + Cm 2 ) ■ 

12 7T 



(175) 



The physical interpretation of this mixing is that the quark condensation has also to be taken into 
account in the Wilsoncocfficicnt of the gluon condensate. The mass of the quarks determines if 
their condensation has to be taken into account in the Wilson coefficient of the gluon condensate. 
A simple rule can be formulated. In the case of heavy quarks the quark condensation has not 
to be taken into account in the Wilson coefficient of the gluon condensate. Moreover, the quark 
condensation can totally be neglected in the corresponding OPE. In the case of light quarks the 
situation is reversed. The quark condensation has to be taken into account in the gluon condensate 
coefficient and in the corresponding OPE. 

These rules can be comprised in the statements that heavy quarks condense mainly trough glu- 
ons while light quarks condense mainly through quarks. The most drastic example is the G 3 
condensate. Light quarks totally decouple from this condensate (see [21] and references herein). 
Decoupling means that the Wilson coefficient is zero. 

The effects just described depend on the relative size of the quark mass compared to the renormal- 
ization point /i, that was introduced to separate long and short distance fluctuations. If the quark 
mass is smaller than /i quark condensation is favored, if it is bigger than /i gluon condensation of 
quarks is favored. 

One way to understand the effect mathematically is to analyze the gluon condensate coefficient. 
The loop integral in the coefficient at large momentum q 2 receives its main contribution from two 
regions of virtual momenta. The first is p 2 ~ q 2 and the second is p 2 ~ m 2 . For small quark masses 
the latter region can't be included in the coefficient Cg since some of the quark propagators soft. 
This piece must be subtracted and absorbed into the matrix element (0|to<7<7|0). While for the 
heavy quarks the propagators are not soft in that region. 



In summary this observation leads to the following expressions. 
1. For pure light quark systems the nonperturbative corrections are 

.physical <0|G 2 |0). 



(176) 



2. For heavy-light systems (to 2 >mi) the corrections are 
G mi (0|migigi |0) + C m2 (0|m 2 g 2g2 |0) + C G ^physical 

(0|G 2 |0) = G roi (0|m 1 g 1 g 1 |0) 
+ (c G + ~^C m ^j (0|G 2 |0) + C m2 ((0\m 2 q 2 q 2 \0) + 1^(0|G 2 |0>) . (177) 

The heavy quark mass expansion for the heavy quark condensate (0 1 hh | 0) = — ^- (0 1 G 2 1 0) + 
... is used to eliminate the C m2 term. The result for the non perturbative corrections is finally 



c,„Js)\nn<h<n\i)) + \ C G + ^~C mi j ;o|G-|o). 



(178) 
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3. For pure heavy quark systems the relevant expression is 



C mi (0|migi<7i|0) + C m2 (0|TO 2 9292|0) + Gg .physical 

(0|G 2 |0)=G G (0|G 2 |0) 
+C mi ^mi<0|g ig i|0) + ~(0|G 2 |0)) + C m2 (<0|m 2 g 2 g 2 |0) + 1^(0|G 2 |0)) . (179) 

Again the heavy quark mass expansion for the heavy quark condensate is used. The result is 

G G (0|G 2 |0). (180) 

It is identical with the one calculated without taking into account the coefficients of the 
miqiQi terms. The conclusion is that the heavy quark condensate has practically no effect 
on the polarization functions. An example is the charmonium system, which in the OPE is 
determined exclusively by gluonic effects. 

In retrospective this section improves the naive picture of the gluon Wilson coefficient. The con- 
densation via quarks or gluons is entangled and can not be treated separately. 



5 The condensates 

In the last sections the Wilson coefficients were in principle determined, but what about the 
condensates? There exists no method to calculate them from first principles. All values found 
in the literature have been extracted from experiments. Some condensates have been directly 
extracted from the data, others involve additional estimations. Hence it is logical to start with 
the ones that have directly been extracted from the data. These are the matrix elements of the 
light quark condensates, some of the four quark condensates, the gluon condensate and the mixed 
gluon condensate. The triple gluon condensate has also been directly extracted from data, but is 
not treated here. 

• (0|tm|0) = -(0.250GelO 3 = -1.5625 • 10" 2 GeV 3 

• (0\dd\0) = -(0.250GeV) 3 = -1.5625 • Kr 2 GeV 3 

• (0|ss|0) = 0.8 x (0|tm|0) = -1.25 • 10" 2 GeV 3 

Which have been computed in the formalism of Gell-Mann, Oaks,Renner [15]. Then there is a set 
of matrix elements which have computed with the help of QSRs. 

• (0|^G^G^|0) = 0.012 GeV A = 1.2 • 10~ 2 GeV 4 

• g(qa^^qG a ^) = M$(qq) = 0.80Ge^ 2 • (-1.5625 • 10- 2 GeF 3 ) = -0.0125GeV 5 

The four quark condensate and the triple gluon condensate have also been computed by using 
QSRs. All other matrix elements have been extracted from the data with less direct methods. In 
the case of the heavy quark condensates the expansion of the condensate reads 



(0\hh\0) = - 




(181) 



The next to leading terms are multiplied by higher powers of Hence, they are suppressed in 
the heavy quark case. The estimates for the heavy quark condensates in first order are: 

• m c = 001.30GeV -> (0|cc|0) = -7.693 • 10~ 4 GeV 3 

• m b = 004.50Gey -> (0|66|0) = -2.2 • 10~ 4 GeV 3 
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• m t = 176.20GeV -> (0|tt|0) = -5.675 • 1CT 6 GeV 3 . 

Thus the heavy quark condensates are negligible compared to the light quark condensates. The 
estimates of higher dimensional quark condensates goes along different lines. The basic idea that 
stands behind the procedure is the assumption of vacuum state dominance. This idea is illustrated 
by factorization of a simplified 4 quark condensate 

(0\qqqq\0) = (0\qqlqq\0) = (0\qq ^ \n)(n\qq\0) = ^(0\qq\n) (n\qq\0) « (0\qq\0) (0\qq\0) . (182) 

n n 

After the insertion of a complete set of states it is assumed that everything except the term with 
the vacuum states can be neglected. Hence, the vacuum is assumed to be the numerically biggest 
contribution. The result is of course only an estimate, but a very handy one. The simplified 
four quark condensate has been factorized in the product of two two-quark condensates. These 
condensates have been estimated as shown before. 

The four quark condensate as defined in (120) has an interior structure which changes the factor- 
ization procedure slightly, but the calculations give insight in an interesting fact and will therefore 
be shown here. 

In many calculations Ti and I^ turn out to be A a 7 M where A a = 2t a and the t a s are the Gell-Mann 
flavor matrices. The condensate for this case is factorized here. Here only quarks of one color are 
considered 

<0|gA 7„OTA o7 „g|0) = (0\q hal fX^q^q^X l qi\0) = ifXh^m^qiMim- (183) 

The interior structure is already factored out, but in this example a non-trivial index structure of 
the quark fields is given and there remains one step before the vacuum saturation assumption can 
be applied. Two kinds of indices have to be considered next. They are the Dirac indices a, /?, r\ 
and v and the flavor indices i ,j,k and /. Vacuum expectation values of operators do not have spin 
or flavor and of course no color. Therefore the flavors and the anti-flavors have to cancel and also 
the operators must be a scalar. To achieve this the fields are grouped in two pairs in each pair is a 
q and q field operator. The indices in these pairs have to be equal in order to cancel the properties 
they correspond to. Hence two combinations are possible. 

• qi, a qj,pqlli 

• Qi.aqi 

Therefore (183) can be written as 

if^rfX'mMqW - MftqjM- (184) 

The double summation in the right summand is cancelled after application of the vacuum states 
and the minus sign stems from the anti-commutator relation that has to be used in order to change 
the position of the quarkfields. As the last step the summation has to be performed. It is achieved 
by employing the relation 



(0\q A q B \0) = 




(185) 



where the quark fields carry a superindcx which collects all possible indices and the quarkfields 
without indices are the summed up fields in which the indices are always equal. The normalization 
factor counts the summands in this sum. In the actual example it is 4 x 3 = 12. Application of 
the relation (185) to (184) results in: 

aa\ii \kk ar/ \ik fi \ki 

% a Jr a <oift„ft„gggzio) - ^ a jr a (oi^tf^io). (186) 
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The corresponding index free notation is 

^2 [* (7m Aa) tr (7 A 'A a ) - tr (7 AJ A Q 7 A1 A a )] (0|«<ra#>. (187) 

At this stage its recommendable to apply the vacuum saturation hypothesis. The result is 

1 16 

[tr ( 7 ^A a ) tr ( 7 ^A a ) - tr (7^7^)] (0|^|0) 2 = -- (0|^|0) 2 (188) 

where also the numerical result is displayed. The full result also has to incorporate the possibility 
of different colors, but the structure of the factorization formula stays the same 

(0|$ri^r 2 ^|0) = \trTxtrT2 - trY 1 T 2 ] (0|^|0) 2 (189) 

as shown in [1]. Here the normalization factor is 4 x 3 x 3 = 36 and (0 1 -0-0 1 0) = (0\(uu + dd + ss)\0) . 
The factorization can also be achieved by using Fierz transformations. An important example for 
a factorization of a four quark condensate is: 

16 16 

a s (q^X a q £ (n 7 ^A a n)) = - a s {qq? = —5- 1-75 • W^GeV 6 = -3.11 • W^GeV 6 (190) 

y y 

n— u,d,s 

(see [15]). Along these lines all quark condensates can be factorized in the approximation where the 
vacuum saturates the whole condensate but if this assumption is valid or not is another question. 
There hadrons in QCD which are very sensitive to the four quark condensate. If the Sum Rules 
are fitted to the data, it turns out that the vacuum saturation hypothesis fails. In order to 
incorporate the violation of the vacuum hypothesis a factor k is introduced. The factorized four 
quark condensate is multiplied by this factor 

(O^r^r^lO) =k~ [trTrfrTi - trT^} (O|0/>|O> 2 . (191) 

This factor takes into account the violation of the vacuum saturation hypothesis. If k = 1 the 
hypothesis holds, if it is not the hypothesis is violated. The four quark condensate is still subject 
to new estimations. The value for k seems to change steadily [16, 36]. However, in this work the 
validity of the vacuum saturation hypothesis is assumed. 

Estimates of higher gluon condensates exist, but are not reliable. Therefore, they have been ex- 
tracted from the data. 

The last remaining condensate on the list is giqcr^ 1 '^qG®„) here also no satisfactory method for 
the estimation or calculation exists and this condensate has again been extracted from data [16]. 
In conclusion it can be said that in contrast to the Wilson coefficients the knowledge of the con- 
densates is poor and relatively high errors are connected with their values. Clearly, the reason for 
this is that they involve non-perturbative physics, while the Wilson coefficients can be computed 
perturbatively. 

6 The dispersion relations 

The central elements in QCD Sum Rules are the OPE and the dispersion relations. Above the OPE 
has been introduced and here the dispersion relations will be reviewed. First the basic formulas 
are introduced and different regularization methods are discussed. It turns out that everything 
is based on only one dispersion relation. The derivation of this relation will be explained in the 
subsequent section. 
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6.1 The formulae used 



The relation on which all other relations are based on is given by 



Imll(s) 



ds 



(192) 



and states that a correlator is known if its imaginary part along the positive real axis is known. In 
the derivation only causality has been assumed. In practical applications (192) has to be checked 
for convergence. ImH(s) can be very steep and the integral can be divergent. In such a case the 
integral has to be subtracted. Therefore, a Taylor expansion of the correlator is performed 



n( q 2 ) =n(o) + 



q< + ... 



+ 



(n-l)\ 



n-1 



it n (« 2 ) 



( q r- 1 + (q2) 



J q 2 =0 



2 ) n f ImU( 

7T J S n (s — 



q 2 ) 



ds. 



(193) 



The first n-1 terms are the expansion while the last term is the remainder term of the expansion. 
Thus, the integral is split in two parts a polynomial in q 2 and a subtracted dispersion integral. 
Divergences of the correlator are eliminated by renormalizing the coefficients of the polynomial. 
For a logarithmic divergence one substitution is sufficient 



n(<z 2 )=n ( o) + ^j 



ImU(s) 



ds. 



s(s — q 2 ) 

After renormalization of the amplitude ensures that II (0) is either finite or 0. If II (0) = 

ImU(s) 



s(s — q 2 ) 



ds. 



(194) 



(195) 



The divergent integral in (192) is finite after one subtraction due to the ^ factor. Another way 
of treating such integrals is the method of power moments M n (g 2 ) , also called Hilbert moments. 
These moments are derivatives of the dispersion relation (192) with respect to q 2 



df U{q)= -JTs-qW dS 



(196) 



d 

~d4 2 



2 \ _ 2 / Imll(s) 



(s - q 2 f 



ds 



(197) 



dq 2 



n(? 2 ) 



2 • 3 f ImIL(s) 



TT J (s-q 2 ) 4 



ds 



(198) 



Or in compact form: 



d Vtw 2\ 2-3-4 f ImU(s) , 

5?) n <«>=— J (H^ ds - (199) 



M -^'h{i)'^=U<^^ d ' n> - 1 - (200 » 
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Figure 8: The analytic structure of the correlators used in this work. On the positive real axis 
above the threshold begins a cut on which additional poles occur. The cut corresponds to physical 
processes. 

At q 2 = M n coincides with the coefficients of the subtracted dispersion relation. This time not 
the correlator itself is computed but its derivatives, this approach will also prove useful in further 
considerations. Again the convergence of the integral is improved with increasing n until the n is 
reached at which it becomes finite. 

A further method to improve the convergence of the dispersion relation (192) is the Borel trans- 
formation. This method is also using derivatives but not a finite number of them but an infinite 
number. The transformation is accomplished by applying the Borel operator 

® = V^Ty { ~ q2)n (dq^J '^ 2 ^°°' 71 Z ^ = M 2 fixed. (201) 

The action of the Operator is to change the kernel of the integral 

BU (q 2 ) = B- r ImUiS ), d S = - j°° e-^ImU( S )ds. (202) 

7T Jo (S - <T) 7T JO 

Hence, after the Borel transformation high momentum components are exponentially suppressed 
for q 2 < 0. For many QSR applications, the Borel transformation is the optimal choice. 

6.2 Derivation of the dispersion relation 

The analytic properties of the correlators together with the residue theorem [23] form the basis 
of the derivation. In a sense the analytic structure of the correlators used in this work is simple. 
They have a cut and poles on the positive real axis above some threshold. The remaining part 
of the complex plane, over which they are defined is free of singularities (see figure 8). Hence, 
according to the residue theorem an integration over a contour which excludes the shaded region 
in figure 8 vanishes. However, a look at the integral in (192) reveals that the integrand has an 
additional pole at s = q 2 . The residue of this pole is the reason why the dispersion integral is 
n (q 2 ) . Instead of a general derivation of the dispersion relation a simple illustrative example is 
discussed. These illustrations mimic the analytic structure of the correlators that are concerned 
here. The particular example chosen is the square root, which has a cut on the negative real axis. 
The analytic structure of the square root is given in figure 9. The only difference to the correlators 
is that there are no poles on the cut of the square root. Everything that will follow from here on 
is simply an application of the residue theorem, which is given by 

n . . 

J2 Res(f, z k ) = —j f(z)dz. (203) 
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Figure 9: The analytic structure of the square root. The cut of the square root is defined to be on 
the negative real axis. 



In this example the integral equation is given by 



The choice of the integration contour is restricted by three requirements 

1. Exclusion of the singularities along the negative real axis. 

2. Inclusion of the pole given by jzj^- 

3. Reduction of the integral to the form given in (192). 

Thus, the integration contour can be chosen as in figure 10. The next goal is to eliminate the 



Figure 10: The integration contour for the dispersion relation before the limit R — > oo is taken. 

integrations along the circles. Therefore, the limits R — > oo and r — > where R is the radius 
of the big circle and r is the radius of the small circle in figure 10 are taken. In many cases the 
integrand vanishes along the circles in this limit. In this example the integrand does not vanish on 
the circle whose radius had been set to infinity, but a simple subtraction eliminates this problem 
(see equation (193)) 




(204) 



k Im [z] 




— ► 

Re [z] 




(205) 
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On the circles the integral is given by 



Re i( J> - z 



e^dz. (206) 



Hence, the integrand vanishes for R — > oo and R — > 0. Therefore only the contour along the 
negative axis survives after the limit is taken. This is a huge step forward but the integration still 
runs over the real and complex part of the square root whereas the one in (192) runs only over the 
imaginary part. 

The solution of this puzzle is given by another feature of the square root. On the integration 
contour that runs above the negative real axis the square root has the value +zy / }z{ on the one 
below its value is — i\f\z\- Such discontinuities also occur in correlators but in general their real 
part is not zero. However, the square root is purely imaginary on both contours. Moreover, the 
the integrations along the negative real axis can be comprised to a single integration because the 
integration on the contour where the square root is — i\/\z\ runs in the opposite direction to the 
one where the root is +iy / jz|'. Therefore the two integrations are equal and can be added. Thus 
the integration takes the form 



fe=i 



Reatf, z k ) = Z f f -r^-rdz = J f ^^dz V~z - Um . (207) 

™ J-oo Z\Z - Zo) 7T J_ oc Z{Z - Zq) 



Only the left hand side has to be modified. It is given by the residues of the integrand in (207), 
but the integrand has only one singularity inside the integration region. Hence, there is only one 
residue. The calculation of this residue is simple and can be performed as follows, the square root 
has to be expanded around z , the residue is the coefficient of the ^-j^-term 

V^= V^+\^(z-zo)-^^{z-zof + ... (208) 

\fz x/Zn 11 11, , 

- +-— --_(*-*„) + ... . (209) 



Z - Z Z- Zq 2 yfzo' 8 

The analysis of the Laurent expansion that has just been derived shows that the residue is given 
by y/zo and the integration (207) finally shows up to be 

^ (210) 

J-oo Z{Z - Z ) 

and for the example of the root the dispersion relation has been valid. Though, there are a few 
remarks that are necessary in order to extend the derivation to correlators. The integration contour 
for correlators is mirrored with respect to the imaginary axes and then the zero point is shifted to 
the threshold value shown in figure 8. This is done because the positive real axis should be excluded 
from the inside of the integration contour. Hence, it has to be justified why the integration along 
two contours can be replaced along the integration over one contour and why the integration over 
real and imaginary part can be replaced with an integration over the imaginary part. 
The solution was already given when the root was concerned but there is a property of correlators 
that has not been explained yet. Along the integrations observed here the following relations hold 

Im [n (q 2 + ie)] = -Im [U (q 2 - ie)] (211) 
Re[U(q 2 + ie)} = Re [U (q 2 - ie)} (212) 
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where q 2 is real and above the threshold. The case of the imaginary part is already known from 
the root example and the result was that the integrations just add up to one integration. From 212 
it follows that the real parts cancel. Thus only the integration over the imaginary part survives. 
As the last step the residues of the correlators have to be calculated but this is a difficult task 
that is not going to be performed here. Instead of another theoretical derivation some practical 
applications are going to be performed. The integral in (210) is now calculated. 
There are two cases that have to be distinguished. 

1. zo does not lie in the region that has been excluded from the integration. 

2. zq lies in the region that has been excluded from the integration. 

The first case is the one which is interesting in applications. Therefore it is going to be explained 
for the root example. Hence in the first case the pole is located inside the integration contour 
indicated in figure (10). Then relation (210) can be used and only a real integration has to be 
performed 

V* = ^ [° -^-,dz = -i^ ( ArTanh L^X) = t Vf° ( = ^ (213 ) 

Thus, the dispersion relation in subtracted form is verified in the case where zq does not lie on the 
cut of the square root. 

The lesson that can be learned from this section is that a correlator is fully determined by its 
imaginary part above some threshold. In the following sections this part of the correlator will be 
identified with the non-perturbative resonance physics, which cannot be calculated by perturbative 
QCD. Hence, if the full imaginary part is known and the correlator is calculated at some point 
using a dispersion relation, the correlator at this point also includes non-perturbative information. 
Dispersion relations are long known and have been used in electrodynamics under the name 
Kramcrs-Kroning relations. A nice review from a modern perspective is given in [25]. 

7 The QCD Sum Rule (QSR) Method 

In the preceding sections the operator product expansion and the dispersion relations were dis- 
cussed. By combining those tools one obtains the so called QCD Sum Rules. The QSR allow 
calculations of properties which are connected to n-point correlators, in a sense these correlators 
can be calculated by QSR. This thesis is restricted to correlators where the external states, are 
vacuum states, but the QSR can deal with arbitrary external states, addressing e.g. the in-medium 
properties of hadrons. 

The objects which are investigated are mainly 2-point correlators. Such correlators contain much 
physical information. The first example for a 2-point correlator that occurs in QFT is the propaga- 
tor, for example (0|T {q(0), q(x)} |0). The propagator has a clear physical interpretation. It is the 
amplitude for the propagation of a particle from the space-time point to x. In the literature [22] 
two cases are analyzed. One without interaction and one with interaction. Without interaction 
everything gets simpler and the propagator is fully calculable and has just one singularity. In 
momentum space it is a pole connected with the mass of the particle. This pole occurs when the 
particle is on shell q 2 = m 2 . With interactions the singularity structure of the propagator is much 
richer. In addition to the pole corresponding to the mass of the particle, poles corresponding to 
bound states and/or resonances occur. 

With QCD Sum Rules one can calculate the properties of bound states by matching the correlation 
functions at a scale between the perturbative and non-perturbative regimes. However, the poles 
can not be calculated in a perturbative approach to the problem. Hence, non-perturbative methods 
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are needed. There is a domain where the OPE agrees with the correlation function. In this domain 
where the correlators have no poles. Thus, an ansatz for the imaginary part of the propagator can 
be made. Furthermore, the space-like region can be connected to the time-like region, where the 
physical singularities of the correlator are located, by means of a dispersion relation. Then the 
propagator on the real axis below threshold can be calculated from the ansatz using a dispersion 
relation. The result can be matched with the OPE in the space like region. This is how QCD Sum 
Rules work. 

Such a procedure allows a comparison between phenomenology and theory because for some cor- 
relators the imaginary parts can be measured. Therefore, in such cases one side of the dispersion 
relation is given by the phenomenology and the other one by theory. These remarks are charac- 
teristic for the way in which QSRs are calculated and evaluated. 



7.1 The connection between the correlator and the spectral density 

There is a direct connection between a spectral function, sometimes called spectral density, and 
the correlator in the corresponding channel. The vector-current correlator is of particular interest. 
It is factorize into a Lorentz tensor and a Lorentz scalar part 

IW«„) = [g^ - U(q 2 ) = (g^q 2 - q^) p(q 2 ). (214) 

The spectral function is proportional to the imaginary part of II (q 2 ). The spectral function in the 
vector channel can be extracted from data on e + e~- or annihilation [22]. In other channels 

no direct measurements exist. 



7.2 The model for the 2-point correlator 

As discussed in the previous section, a QCD Sum Rule is a dispersion relation connecting the 
phcnomcnological spectral density with the OPE 

i r> im[u ph ( S )] ds = QpE 

7T Jo (s-q 2 ) 

The momentum domain has to be restricted to large space-like momenta. Hence, the substitution 
Q 2 = —q 2 is convenient to avoid the unimportant minus sign. 

Obviously a QSR can be used in two ways. The first one is to define the right hand side as given 
and to calculate the properties of the left hand side or vice versa. Both possibilities are used 
and have to be used, as will be discussed below. The evaluation of a QSR is a chapter for itself 
and is best illustrated through examples. The classical way to use QSR is to use an ansatz for 
the phcnomcnological part given by a narrow resonance and a continuum, where the resonance is 
given by a delta function and the continuum by a step function 

ttpheno(s) = (0|j|^) 2 <5(s - M res ) + const. ■ 6{s - t c ). (216) 

If a tiny width is given to the <5-functions the phenomcnological part in the form sketched in figure 
11. Obviously such a model is very simplified, since the particle that is represented by the 6- 
resonance is usually a resonance with a non-zero width. Thus, a resonance with a finite lifetime is 
represented by a stable particle. Moreover the radial excitations are all missing. The step function 
is an approximation for the continuum. The first step to make calculations handy is to use the 
imaginary part of the perturbative expression for the correlator in the domain where q 2 is positive 
instead of a step function. Especially because the perturbative expression is a part of the OPE 
which belongs to the system and is therefore known. Others improvements will be shown later. 



50 



ri(s) 















s 



Figure 11: Model for the phenomenological part of the sum rules. 



The right hand side is the theoretical side which is given by the OPE. However, the OPE can not 
be calculated exactly, it has to be truncated. In Section 3.4 the criteria for convergence of the OPE 
were discussed. Furthermore, the terms in the truncated OPE are not known exactly, but to a 
certain order in as or to a certain accuracy of estimation, as discussed sections 4 and 5. Finally the 
exact Wilsonian OPE is not employed but rather an approximation, the so called practical OPE or 
SVZ expansion (see section 3.7). Hence, also the theoretical side of the QSR contains uncertainties 
and simplifications. However, QSRs have been and are heavily used in non-perturbative physics 
with big success and with a satisfactory accuracy. In many cases QSR arc the only tools available 
on the market. 

A well known applications of QSRs are the p- and the J/ip-meson. The analysis of these two 
systems is instructive and will be presented in the following section. 



8 The J/ip QSR as a classical QSR 

In this section QSRs for the J/V'-meson are employed to determine the gluon condensate as(G\^ v G^\0) 
8.1 The spectrum of cc-mesons 
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Table 4: The known charmonium spectrum. The excited states correspond to radial excitations. 



In table 4 the known spectrum of cc-meson (Charmonium) is shown. In the vector channel 
(J PC = 1 ) the ground state and several radial excitations are known while in the other chan- 
nels only the ground states are known. Such a situation is characteristic for hadron spectroscopy. 
In many cases only the ground state corresponding to given quantum numbers is known. Therefore 
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the calculation of the ground states properties would already be a remarkable success and for this 
system QSR do the job very well. SVZ realized during the work to their paper [3] that cc-mesons 
are bound exclusively by gluons. Thus, an attempt to calculate the mass of a cc-meson with a 
QSR in which an OPE is used which contains only the perturbative and the lowest order gluon 
condensate term is reasonable, although the rules of section 3.4 say that all terms up to dimension 
six should be included. It turns out that the attempt is successfull, but to achieve the goal a 
sacrifice has to be made. Section 5 states the value of the gluon condensate (0| ^-G° G° \Q) and 
as the source of it the J sum rule. The gluon condensate has not been calculated up to now 
and therefore it has to be measured. The cc-mesons are the best systems for it because of the 
structure of the OPE and because the spectrum of the J/'i'-meson can and has been measured to 
a high accuracy. 



8.1.1 The phenomenological part of the sum rules 

The spectral function of the J/'i'-meson and its radial excitations has been measured in 
annihilation (see for example [22]). The resonances can be approximated by 5- functions and the 
spectral function shown in figure 12. The parameterization of the spectral function is very intuitive: 



n(s) 



Figure 12: Sketch of the spectral function of the J/'i'-meson as it is known today. 



ImU(s)= Yl (V\j^res) 2 5( S -ml) + —h+^^)e(s-t c ). (217) 

resonances ^ ' 

In the sum rule for the ground state of the states with J p — 1~ an approximation to the spectrum 
is used where all excited states are neglected and the ground state is approximated by a <5-function. 
Thus, the effective spectral function is given by figure 13. The evaluation of the sum rule will show 
that this approximation is reasonable. In the evaluation process the contribution of the neglected 
parts of the sum rules are suppressed. 

8.1.2 The theoretical (QCD) part of the sum rules 

The current which approximates the J/'i'-meson is given by 

j^x)=c{x) ltl c{x) (218) 
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ri(s) 



Figure 13: Approximation for the spectral function of the J/^-meson used in the sum rules for 
the calculation of the gluon condensate. 



and has exactly the quark structure and quantum numbers of the J/ 4'. In the OPE for the 2-point 
correlator of this current only two terms are kept 

i J dtxe-V^faWMO)) = (g^q 2 - q»q v ) [c pert ■ I + C& ■ (OI^G^JO}] , (219) 

where the coefficients are given by: 
1. pcrturbative contribution C pert 

The approximation to the perturbative coefficient used here contains the bare loop and as 
contributions 



Im 




■ + ■ 





= Im [C pert (s)} = Q{s - 4m 2 )-"(«)(3 - v 2 {s)) 

87T 



v(s)+3 



2v(s) 



3 

47T 



v(s) = 



Am 2 
— • 



(220) 



Dispersion relations can be used to calculate the real part of the amplitude. It is instructive 
to plot the perturbative part of the OPE, the plot is given in figure 14 and shows only 
the amplitude corresponding to the first diagram in (220), the one loop diagram. Higher 
order corrections do not change the overall structure of the plot. Thus, the lowest order 
contribution shows the important features of the perturbative part of the 2-point correlator. 
The imaginary part is going to be used as an approximation for the continuum in the model 
for the spectral function. The amplitude is real below the threshold, above threshold the 
imaginary part is non-zero. 
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Figure 14: The one loop contribution as a function of the total momentum squared. The dashed 
line is the real part and the full line the imaginary part. 



2. nonperturbative contribution Cqi 

The only nonperturbative contribution stems from the gluon condensate term. The coef- 
ficient is determined by the bare loop with two external gluon legs (see section 4) 




In contrast to the diagrammatic description the coefficient is not a simple diagram. Its plot 
is given in figure 15. The coefficient Cc(q) in (221) is valid in the limit \q\ ^ 0. Hence, it 
is expected to find an unreasonable behavior of the coefficient if q approaches zero. In figure 
16 the coefficient is plotted in the momentum region close to zero. There a highly oscillating 
behavior is found, which signals the invalidity of the coefficient in that momentum domain. 
Unfortunately, this effect is not visible in figure 15. 

Thus, the sum rule has been constructed and can now be evaluated. Momentum Sum Rules are 
employed for this task, see section 6 for the basis that is needed for further considerations. 
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Figure 15: Plot ol Cq (q 2 )? the dashed line is the real part and the full line the imaginary part. 
The singularities above q 2 — occur at 4m 2 . 
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Figure 16: Plot of Cq (q 2 ) showing highly oscillating singularities at q 2 = 0, the dashed line is the 
real part and the full line the imaginary part. 
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8.1.3 Evaluation of the Sum Rule 

If (217) is used in the dispersion the following equation is obtained 

1 ( d Vn(c 2 )-V mti^i i 1 f 1 + Qis t)ds n> 1 (222) 

n\ \ dQ 2 ) 11 [Q >~^ s (m 2 fl + Q2)(n+D + 47 r 2 J ( s + Q 2 )(«+D ° (S tc)dS " " 1 (222 > 

where n has to be bigger than to make the integral convergent. This relation allows the calculation 
of the J/ip mass. The mass of the ground state with the quantum numbers J p = 1~. A simple 
conversion will show that this is true 

- 5 + i / (^f^ e(s - (223) 

The M n 's are called the Moments. First of all the term corresponding to the J/i/j-resonance is 
factored out 

M n (Q 2 ) = W>'/*> a 



K /5/ , + Q 2 ) ( " +1) 



The quantity 



fm 2 , + 2 ) {n+1) 

d ™ W ) = JKTa l„ \2 



(0|.7 M |nj/v.) i 



(0|j„k„> 2 , 1 ' 



res>J/ip 

is convenient because £ n is suppressed for big n 



0, n — ► 00 



m re S +Q 2 J 

n+l 

) — >0, n — > 00, s>t c . (226) 



This means that: 



Now the ratio 



s + Q 2 



M n _! m 2 ^ + g 2 1 + 5, 



-1 



is considered. Depending on the channel, the factor 1 1 ^ s Sn - can sometimes be replaced by 1. The 
spectrum of J p = 1~ cc- meson is such a case. This can be tested because the spectrum has been 
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measured. This is equivalent to the replacement of the spectrum in figure 12 with the one of 
figure 13. If the r„'s are also computed by the OPE, the right hand side of the dispersion relation, 
theoretical predictions are possible. 

A simple example would be the calculation of the J/ifi mass. The c-quark mass is known sufficiently 
well to expose bounds on it and therefore there is only one really free parameter remaining, the 
gluon condensate. Hence, the mass calculation is possible, but the c-quark mass and the gluon 
condensate have to be fitted. Thus, the gluon condensate can be calculated by adjusting it to a value 
which reproduces the mass of the J/i/j. The result is shown in figure 17. At small n the continuum 

4.5 t j 

4 - 
3.5- 

3 ~ ' " ■ s ■ i~~ 

raj/f in GeV 2.5- , - 

2 - ; 

1.5- \ 
1 - 

2 4 6 8 10 12 
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Figure 17: Plot of the moments M n (~lGeV 2 ) belonging to the QSR of the J/vp-meson. The full 
line is the measured mass of the J/^S, while the dots are the moments at a squared momentum of -1. 
The zeroth and the first momentum do not converge due to a dispersion integral used to calculate 
the perturbative part of the Wilson coefficient. In the domain of the plateau, the contributions of 
the perturbative and non-perturbative terms in the OPE are optimized. 

is not suppressed strongly enough and gives the dominant contribution to the phenomenological 
side. On the theoretical side the perturbative part of the OPE gives the main contribution to 
the Moments. Fortunately, with growing n the continuum is suppressed and the resonance gives 
the dominant contribution to the phenomenological side. Unfortunately, with further increasing 
n the contribution of the non-perturbative terms to the theoretical side becomes the dominant 
contribution. The plateau lies in the region where the resonance dominates the phenomenological 
part of the sum rules and the perturbative part of the OPE dominates the theoretical side. 
The plateau determines the mass of the J/ip. In the n region of the plateau, the contribution of the 
perturbative term to the moments dominates and the non-perturbative term of the OPE (219) is 
smaller by factors of 10. The plateau starts when the non-perturbative contribution to the moments 
is a 100 times smaller than the perturbative and ends when the non-perturbative contribution is 
equal to the perturbative one. Therefore left to the plateau everything is perturbation theory and 
there are no bound states. Hence, there is no plateau. While on the right of the plateau the 
non-perturbative terms dominate, but they are not known accurately enough and therefore the 
sum rule does not converge any more. In the domain of the plateau the convergence of the sum 
rule is optimized between convergence of perturbation theory and the lack of knowledge concerning 
the higher terms in the OPE. The c-quark mass and the gluon condensate are fitted in order to 
reproduce the mass of the system. The result is m c = 1.4 GeV and the gluon condensate has the 
value quoted in 5, in this example the moments have been evaluated at q 2 = — 1 GeV 2 other values 
are of course also allowed. 
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Historically this was the starting point for the calculations of QSR [1]. The gluon condensate 
has up to now never been calculated from first principles. It is multiplied by the strong coupling 
constant in order to form a renormalization group invariant. 

Now this value of the gluon condensate can be used to calculate QSR for other systems. Reinders, 
Rubinstein and Yazaki performed a rigorous test of the gluon condensate on the cc spectrum. 
They repeated the calculation just performed for every quantum number J p for which the ground 
state was measured. These calculations confirmed the gluon condensate given in section 5 (see [5]). 
Moreover many other systems have been analyzed during these times. 

Many systems containing two quarks with nearly equal masses have been calculated often using 
Borel transformed QSR (Borel Sum Rules) in order to improve the convergence of the Sum Rules. 
The mass of the p-meson given by QSR is remarkably close to the value measured and the plateau 
is very long, (see [1]). In the period from 1978 to 1988 hadrons seemed to be describable in the 
way just outlined, but with proceeding time problems occurred. They occurred in exotic systems 
like hybrids or glue balls and in heavy-light systems. In the case of the D-meson, the sum rules 
converged very badly, and a broad plateau is not given. 

9 Borel transformation 

Section 8 discussed the general pattern how QSR are used. It is allways a dispersion relation 
consisting of a phenomenological and a theoretical (QCD) part. These parts are the left and the 
right hand side of the dispersion relation (see for example (215)). The phenomenological part of 
the QSR can be determined by the intuition of the user, by input from alternative theoretical 
approaches or from measurements, while the theoretical part is given by the OPE. The input from 
both sources is afflicted with errors. The spectral functions are in most cases approximations 
and can not be directly measured with the exception of spectral functions of the vector-current 
correlators. Furthermore the OPE is also an approximation because it is a truncated expansion. 
This is a truncation on several levels. First, it is an expansion in -i?. Second, the Wilson coefficients 
are approximated. Finally, the condensates are approximate. It is difficult to judge the combined 
effect of the approximations. Nevertheless, the QSR are in many cases very successful. 
The circumstance that ensures the reliability of the QSR is the existence of a working window. If 
the OPE would be known completely for the given spectrum the dots in figure 17 would not change 
much for small n but they would drastically change for the n to the right of the plateau. They 
would stay close to the line which shows the value of the mass of the system which is observed 
and would approach it. There exist exactly solvable models which substantiate this statement, 
one is given in [15]. In the ideal case just mentioned the OPE would not be truncated in any 
instance and its the truncation that is responsible for the break down of the plateau. To the right 
of the plateau the non-perturbative terms in the OPE are more important than the perturbative 
ones, but they are not known well enough for reliable calculations and therefore the results become 
unreliable and the plateau breaks down. The working window is therefore defined by the criteria 
of the dominance of perturbation theory or in other words the perturbative part of the OPE. This 
part is known to a sufficient accuracy and the errors that stem from the non-perturbative part are 
negligible as long as the perturbative part dominates in the evaluation of the sum rule. A common 
abbreviation used for the requirement is that asymptotic freedom still holds in the domain where 
QSR are used, which can be very misleading. The existence of such a windows is a hypothesis and 
has to be checked from case to case. 

Many publications just use a guess for the spectral function. Hence, the situation can arise that 
the OPE is under control, but the form of the spectrum is incorrect. In this scenario the sum 
rules may be completely arbitrary. There may be a plateau, but the properties corresponding to 
it are incorrect or the plateau can simply be missing and no statement is made even if the OPE is 
correct. In this case there would be no working window and the OPE would be discarded. 
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These two cases have always to be concerned in the evaluation of a sum rule. Hence, a possibility 
to reduce the sensitivity of the sum rules would be more than welcome, and sometimes mandatory. 
Fortunately such a possibility exists, the Borel transformation which was introduced in section 6 
in order to improve the convergence of a dispersion integral. 
The Borel transformation improves the QSR in three ways. 

1. The dispersion integral is regularized because the Borel transformed sum rule is exponentially 
suppressed at high energies. 

2. The exponential suppression also minimizes the influence of excited states. 

3. The unknown terms in the OPE are suppressed. 

Hence, the Borel transformation improves the convergence of the sum rule if the ansatz for the 
spectral function is a good approximation for the corresponding OPE. The Borel transformation 
is performed by applying the operator 

B = . 1 ,A Q 2 ) n , Q 2 ^oo, n^oo, — = M 2 fixed (229) 



(n-l)P^ ' \ dQ 2 J ' ^ ' ' n 

where M 2 is the so called Borel parameter. The improving features of the OPE can be illustrated 
by two simple calculations. 

1. Application of B to ^ — ^ 5 



s + Q 2 J (n-iy.^ ' V dQ 2 J s + Q 2 (n-l)! v ^ 7 (s + Q 2 )«+! 
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v n ' 



e~^. (230) 



An exponential kernel in the dispersion integral suppresses the integrand for large s and 
makes the integral finite. At the same time the spectral function is suppressed for large s. 



2. Application of B to 
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Clearly the Borel transformation leads to an improved convergence of the OPE. 
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For the perturbative part of the OPE the Borel transformation of the logarithm is needed- 
Application of B to In (Q 2 ) yields 

8 ~ MO - j^WV (-^)'W> - ^(Qr (-^) - -i. (232) 

From the point of view of mathematical consistency it would be better if the evaluation of a QSR 
would work by moment sum rules where the improving features are much weaker. Generally a sum 
rule is more reliable if the moment sum rule converges. However, there are sum rules where the 
moment sum rules do not converge but the Borel sum rules do. A nice pedagogical review to this 
topic is given in [2] , where a sum rule is constructed which does not converge without being Borel 
transformed. 



10 Heavy- Light Systems 



In the sector of heavy-light systems great progress has been made in recent times. A major break 
through was a discovery made in 2004. In the sector of particels with charm and strangeness two 
new particels where found, the D*(2317) ± with quantum numbers J p — + and the £) s (2460) ± 
with quantum numbers J p = 1 + . This discovery confirmed theoretical predictions already made 
in 1992 [27]. Together with the already known states Df with J p = 0" and with J p = 1" a 
hypermultiplett is formed which consist of two doublets (see figure 18). The doublets are given by 



D*±, J p = 1- 
J 



£) s (2460) ± , J p = 1+ 



P 



D±,J p = 0- 



4) D*(2317) ± , J p = 0+ 



Figure 18: Hypermultiplett build up by the most prominent _D s -states. 



two points which are connected via a dashed line, the coordinates of the dots are their quantum 
numbers J p . The quantum numbers are determined by the total angular momentum, the orbital 
angular momentum and the coupling of the orbital angular momentum and the spins. The parity of 
a D-meson is given by P = (— 1) , where I is the quantum number of orbital angular momentum 
of the valence quarks. Hence negative parity is given by even I and positive by odd I. Therefore 
the orbital angular momentum of the left doublet is odd and of the right one even. Moreover the 
particels of figure 18 should be as lowest lying in the spectrum as possible. Hence, the smallest 
I possible is assumed for them. This results in I = for the left doublet and I — 1 for the right 
doublet. The spin configuration for the left doublet is antiparallel for the the 0~ state and parallel 
for the 1~ state, while for the right doublet the spins are parallel for the + state and anti-parallel 
for the 1 + state. 

A descriptive reason for the existence of such particles is given by symmetry arguments. There are 
two symmetries that determine the structure of the hypermultiplett, the chiral symmetry and the 
heavy quark symmetry. In the limit where the symmetries are exact and not spontaneously broken 
the four particels of figure 18 are degenerate in mass. Both symmetries are explicitly broken, and 
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restored in two totally different limits. Chiral symmetry in the limit of vanishing mass and heavy 
quark symmetry in the limit of infinitely large mass. If chiral symmetry is assumed to be exact 
but spontaneously broken for the light quarks and heavy quark symmetry to be exact for the 
heavy quarks, the main contributions to the mass degeneracy and therefore the overall structure 
of the masses of the hypermultiplet can be interpreted in terms of the light quark condensates, as 
discussed by [28]. 

Historically the states with negative parity have been known for a long time before the states with 
positive parity were discovered. A theoretical link to the existence of the positive parity states is 
given by chiral symmetry arguments. The Lagrangian which describes only the negative parity 
states of a heavy-light system does not posses chiral symmetry in the limit of vanishing light quark 
masses. However, a Lagrangian which describes the hypermultiplctt where both, the negative and 
the positive parity states are included, exhibits chiral symmetry in the chiral limit. Such symmetry 
arguments do not prove the existence of the positive parity states, but they are strong arguments 
and show that theory is on the right way. 

Heavy quark symmetry links the J = states directly to the J = 1 states. In the heavy quark 
effective theory, the theory in the limit of infinitely heavy quarks, the spin decouples from the 
dynamics of the system. In that limit a heavy-light system with parallel spins is completely 
degenerate with a heavy-light system with anti-parallel spins. Hence, the states in the doublets 
are linked to each other. 

These considerations have raised a lot of interest for the spectroscopy of £> s -mesons. One expects 
similar arguments to hold also for other heavy-light systems. However, the spectrum can not be 
reproduced exactly by such methods, because only the main contributions to the mass spectrum 
are included and the fine structure of the spectrum is missing. One method to calculate the mass 
spectrum more accurately is provided by the QSR. In the following, the QSR in the heavy quark 
effective theory (HQET) and the results for the masses of the particels in the hypermultiplctt 
are shown. After that the full relativistic theory is going to be reviewed as it has been used 
and accepted for a long time. The shortcomings of such approaches concerning the calculation of 
hadron properties will lead to the next topic. 



10.1 Heavy flavor sum rules for the D-meson 

All of the arguments just made can be checked by QSR calculations. This section concerns the 
question whether the doublets are degenerate in the heavy-quark limit or not. The masses of the 
J p = + and J p = 1+ and those of the J p = 0~ and J p = 1~ states respectively are expected to 
be degenerate in the heavy quark limit. 

In the remainder of the section a mini review of the HQET is given in order to illustrate the most 
important computation tools. The basis of all calculations is the Lagrangian of the system which 
is in the actual case split into two parts. The light quark dynamics is described by the ordinary 
QCD Lagrangian while the heavy quark dynamics is described the heavy quark effective theory 
Lagrangian 

-Sf = J^QCD + ^HQET = flight + ^heavy (233) 

Those Lagrangians are given by 

Xqcd =q{0- m q ) g - \Gl v G% ^ H qet = hilph = hiv^h. (234) 

All that is necessary for further considerations are the Feynman rules of the theory. The propagator 
for the heavy quark is given by 

-L-i±± = iU + l) (235) 
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and the vertex for the heavy quark by 



ig s v^t a . (236) 

The rules for the light quark are the ones of ordinary QCD. The heavy quark vertex does not 
contain any Dirac matrices. Hence the interaction is independent of the spin. The momentum in 
the heavy quark theory is split up in two parts 

9 M = fUheayyVf, + Pf, = 1, (237) 

where mv a is the on-shell part of the momentum and p^ the off-shell part, p^ is also called 
the residual momentum. The splitting of the heavy quark momentum is useful because the heavy 
quark is expected to be non-relativistic. The huge mass is much larger than the typical momentum 
in a hadron. Hence, higher order terms in p 11 can be neglected. The dependence on t> M and p^ can 
be comprised in a dependence on uj = 2v ll p tl . 

The mass of the heavy quark does not appear in the Lagrangian of the HQET, but in the finite 
quark mass corrections. A reflection of the "missing" mass parameter is found on the hadron 
side of the sum rule. The spectral function runs over w and the resonances occur at the points 
2A = Mhi-system ~ mheavy quark, where 2A is the value where the heavy quark mass is set to 
infinity. In this limit the heavy light system also contains a infinitely heavy quark. 
In a effective theory the fields get modified which is reflected by the Lagrangian (234) where the 
heavy quark fields h occur, but also the current operators change. Unfortunately the situation here 
gets very difficult in comparison with the situation in QCD. The currents are not unique in heavy 
flavor sum rules (see for example [29]). However only the following currents in table 5 are used in 
the sections below. These currents are used to compute sum rules for the D-meson hypermultiplctt 





current 


0+ 


qh 


o- 


qihh 


1+ 


975 (7m -i>Vp)h 


1- 





Table 5: Table with the currents which are used for the approximation of the analyzed mesons in 
the HQET. 

which reproduce the expected mass degeneracy in the heavy quark limit. The calculations in the 
HQET are illustrated in the following. 

10.1.1 Calculation of the heavy-light quark loop 

The calculation for the pseudo-scalar and scalar current correlator are demonstrated. In the scalar 
case the current is 

j{x) = q(x)h Q {x), j\x) = h Q {x)q(x). (238) 
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The current correlator is given by 

n p ert(q 2 ) = i J d 4 xe^(0\T(j(x),f(x)) \0)=ij d 4 xe i " x {0\q(x)h Q (x)h Q ( x)q(x)\0) 
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J (2tt)4 



-tr 
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J (2^) 4 
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r) 4 (p2-m 2 )2(t;M(g M +p M )) V (2^) 4 (p 2 - m 2 ) (w + 2^) ' (239) 

where the double line represents the propagator in the heavy quark effective theory. The expression 
for the correlator is renormalized by using the MS'-scheme. Therefore it has to be regularized by 
dimensional regularization. To begin a Feynman parameter is introduced and the integration is 
transformed to a spheric symmetric integration 
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(240) 



with 1^ = p^ + v^^^T and A = (^4r) 2 + m2 ~ ^r w - ^he ^ integral is extended to d dimensions 
and the integration is performed 
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After the subtraction and in the limit d — > 4 the following expression is derived 



Ilpertiq 2 ) = -7-2 
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(242) 



This form of the amplitude is very inconvenient because of the Feynman parameter integral. One 
method to get rid of it is to consider only the imaginary part of the amplitude. Which is sufficient 
for the derivation of the whole amplitude (see section 6) . The logarithm determines the imaginary 



part of the correlator. It is imaginary if the argument is negative 



< 0. This domain is located 



in the interval between the roots of A(q 2 ). Hence, the imaginary part of the correlator is 
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(244) 
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Consequently the imaginary part of the correlator is determined, but the color structure of the 
problem has not been considered. A factor 3 has to be added 

o rx 2 1-x I m 2 

Imn pert (q 2 ) = — / dx ^ = —(uj- 2m)VoJ 2 -4m 2 . (245) 

47T J X1 X Z 87T 

The next step is to calculate the two point correlator of the pseudo-scalar current correlator, where 
the current is 

j(x) = q(x)i^ 5 h Q (x), j ] (x) = -h Q (x)ij 5 q(x) (246) 




(2tt) 4 (p 2 - to 2 ) (w + 2v^Pfj) ' 



The only difference between (239) and (247) is the minus sign in front of v^p^. The result for 
n pert ,5(<7 2 ) is 

ImU pert ^(q 2 ) = ^-(lj + 2m)Vw 2 -4m 2 . (248) 

This is an example for a interesting rule. The results for chiral partners are obtained from each 
other by exchanging to with — m. 

10.1.2 Calculation of the quark condensate coefficient Cq q 

Again the example of the pseudo-scalar current is considered (see table 5) 



/i 1 _ , 1 

d 4 xe lqx (j), s\q(x)h Q (x)h Q (x)q(x)\p,s) = 




= iU { p,s)^^t±±u M . (249) 

The results of section 4.1 are used to derive 

Cq q (qq) = -— ( u , 1 , r - P>i + m ) s )u(p, s) = -— - P ^ + m u(p, s)u(p, s). (250) 

The expression 



1 / w^p^, + m \ 
to \w + J 



(251) 
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has to be averaged over the fourdimensional Euclidean angle as it is shown in section 4.1. The 
result is in the limit m — > is 

Cqq = --. (252) 

The scalar current correlator is determined by the expression 



Cqg{qq) = 




(253) 



q 4 V 



The change in the expression in comparison with the pseudo-scalar case is just the minus sign in 
front of the m-term. The result in the limit m — > is 

Cqq = ~ (254) 

where the difference to the pseudo-scalar case (252) is just a flip in the sign of the coefficient. 



10.1.3 The Borel transformed heavy flavor sum rule 

The Borel transformed heavy flavor sum rule for the D — mesons with isospin or strangeness, in 
the limit mu g ht = reads for the positive parity doublet 

(0|i|0) 2 e- 2 ? = ^ £ u?e-*du, + »^ - ^M 2 (0k-<z|0) (255) 

and for the negative parity doublet 

(0|i|0) 2 e- 2 ? = JL £' Je-tfa - Em + _l_ Mo 2 (0|^|0). (256) 

T is the Borel parameter. The currents employed for the calculation of the sum rule have an 
additional factor (see [29]). Hence, in the heavy quark limit the masses of the doublets are 
degenerated. 

The result of the sum rules evaluation in the channel with isospin is A = (l.OO)GeV for the positive 
parity doublet and A = (0.95)GeV for the negative parity doublet. Hence, the positive parity states 
are supposed to be heavier than the negative parity states. Thus, the prediction of the mass for 
the D — mesons with positive parity is mo ~ A + m c = 2.3 GeV . In the case of the negative 
parity states mo ~ A + m c = 2.25 GeV. The Borel curves are given in figure 19, the curves for 
the positive and negative parity states stay nearly equidistant for T > 1 . Dai et al. where the first 
ones to analyze such sum rules (255), but they just analyzed the positive parity states. However, 
their value for the threshold oj c was used for the calculations shown above. In their publication 
they obtained A = 1.05GeV with uo c — 2.65GeV. The difference between the calculations above 
and their calculations is probably caused by a difference in the condensates. Unfortunately they 
do not give the values they used for the condensates. 

However, in this calculation the values for the quark condensate, as given in section 5, has been 
used. Hence, in the negative parity channel the results can be compared with the masses of the 



65 




Figure 19: This figure shows the Borel curves for A, the mass parameter in the HQET, for the 
states in the hypermultiplet of figure 18. The full line is the curve for the P = + states while the 
dashed line is the curve for the P = — states. 



£>°(1869),L> ± (1864),£>*(2007)° and L>*(2010) ± particles. In the positive parity channel the data 
of the Z? 1 (2420)° and Di(2430)° particles can be compared with the results. Obviously the positive 
parity states lie much closer to the results than the negative parity states. In the P = + case the 
error is about 0.150 GeV while it is about 0.400 GeV. This may be due to the OPE, an improved 
OPE may give better results. 

10.1.4 The mass formula for heavy flavor sum rules 

The intradoublet mass splitting is another topic. It stems from the breaking of the heavy quark 
symmetry. This violation is in QFT formulated as — - — corrections. The Lagrangian with the 
first corrections is given by 



O 



1 



2m Q 2mQ 

where K, is the nonrelativistic kinetic energy operator, defined as 

JC — h v (iDj_) 2 h v 
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'zUf is the first coefficient of the f3 function. 



Taking into account the corrections in the Lagrangian, the meson mass formula in HQET is 
expressed as 



M = m Q + A - 



2m ( , 



(Ai + d m X 2 ) 



(260) 
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where the two additional parameters Ai and A2 at the -7— order are defined by two matrix elements 
2MXi = (M\!C\M), 2d M C mag MX 2 = (M\S\M) (261) 

the constant cLm is spin-related <1m = djj t and d- ( _ 1 - ( = 2ji + 2, dj l + i j t = ■ 
The evaluation of Ai and A2 needs the consideration of the three-point correlator 

S(w,w')=« 2 y d i xd i yexp(ikx - ik'y)(0\Tj\x)O(0)J(y)\0) (262) 

where the operator O can be JC or S and J is still the generic interpolating current. The derivation 
of the formulas can be found in [30] and references herein. 

10.1.5 Shortcomings of the heavy flavor sum rules for the D-mesons 

The preceding sections have shown the characteristics of heavy flavor sum rules and their short- 
comings. The mass formula shows that a lot of additional work has to be done in order to get the 
corrections to the heavy quark limit. Moreover, the operators in the HQET have to be renormal- 
ized in order to get the right results (see [31]). This is also not necessary in the full theory. In 
summary the heavy flavor sum rules can be excluded as the right tool for the calculation of the 
properties of the spectral function, because except for additional work they supply nothing new 
nor an improvement against the sum rules in the full theory. 

Even if the first corrections are taken into account the effective theory does not produce a satisfying 
accuracy (see for example [32]). This is a second argument against heavy flavour sum rules for the 
calculation of properties of spectral functions. 

10.2 The QSR for heavy light-mesons in full QCD 

The first step is to set the boundary conditions for the OPE. In section 3.4 a criteria was formulated 
which allows a classification of the terms in the OPE. Important are all terms containing opera- 
tors of dimension smaller or equal to the dimension of the operator product under consideration, 
d(Oj) < d(jj). All remaining terms are less important. The currents for the heavy-light mesons 
which are analyzed here are given by 

js = (m 2 + mi) q 2 qi = <9 M jy 
jp = (m 2 - mi) q 2 ij 5 qi = d^f A 

iv = qii^qi 

3 A = 52757^1 ■ (263) 

The scalar and pseudo scalar currents have dimension four while the vector and axialvector current 
have dimension three. Hence, all important operators in the scalar or pscudoscalar case have 
dimension less than or equal to eight, while in the vector or axialvector case the important operators 
have dimension of less than or equal to six. In the ideal case the OPE for the particels approximated 
by these currents would contain all of these terms. However, the situation in reality is far from 
that goal. 

Historically the pioneers in the area of mass spectroscopy of heavy light systems with QSR were 
the authors Novikov, Voloshin, Shifman, Vainshtein and Zakharov (SVZ) in the late seventies [4]. 
Moreover the authors Reinders, Rubinstein and Yazaki (RRY) continued their work in the early 
eighties. They computed the first QSR for this problem (see [6] and references therein). They 
truncated the OPE on a very crude level, keeping only the quark and gluon condensates. However 
their work showed that the sum rules for such particels work and they could derive first estimates of 
properties for the B-meson family with moment sum rules. Properties for the D-meson family were 
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first computed by Alicv in 1983. He used an OPE which contained more operators than RRY. This 
lead to predictions of lcptonic decay constants for D-mesons (see [12]). He included, in addition 
to the quark and gluon condensate, the mixed and four quark condensates in the OPE. Narison 
used an OPE which was similar to the one of Aliev for further investigations of heavy-light systems 
in the late eighties [13]. His work showed that it is possible to treat B- and D-mesons with such 
sum rules, but for D-mesons the use of Borel sum rules was necessary. Moment sum rules do not 
converge for the D-meson. Despite the differences in the approximations made by the groups they 
also have things in common. The continuum contribution ca not be neglected as it has been done 
for the charmonium system. This is nothing particular to heavy-light systems. The continuum 
is important in many cases. The p-meson sum rule is a classical example where the continuum 
has to be taken into account. During the nineties the heavy quark effective theory entered the 
realm of QSR but as shown in section 10.1 this approach has not brought much progress to mass 
spectroscopy with QSRs. Thus for many years Narison was the only author who pushed the field 
forward. He always used the sum rules close to the one constructed in [12]. Aliev's sum rule was for 
D-mesons with J p = 0~. Thus all analyzes based on this sum rule where done for 0~ or + states, 
the sum rules for the two particels are connected by a transformation. There have been many 
papers about QSR and heavy-light systems which simply used this sum rule in order to perform 
calculations, but they only applied Narison's work. The situation changed after the discovery of 
the positive parity doublet in the D s -meson family. The interpretation of the nature of such states 
is still open and has been attacked in papers like [18], [20], [33] and [34]. Where Narison's [18] result 
implies that the dominant configuration of the Ds (2317) is a quark anti-quark state, Hayashigaki 
[20] finds, using an unconventional approach that a more complicated state must dominate. In 
both publications two quark currents and the sum rule of the Aliev type are used for the J = 
particels. An important difference is that Hayashigaki uses a sum rule for the J — 1 states, while 
Narison uses HQET arguments to derive the masses of the J — 1 particels from the J = masses. 
Hence, Hayashigaki's calculation is superior in the J = 1 channel. In the remaining publications 
a four quark current was used as the interpolating current for the Ds(2317). However, there still 
remains much work to be done. 

A short review of the sum rules used by Narison and Aliev exhibit the shortcomings of those sum 
rules. The OPE for the sum rules is of course truncated and the Wilson coefficients are calculated 
to a certain order in the strong coupling constant as- Moreover the limit of massless quarks is 
used. A curious point is that various authors do not agree on the form of the OPE. Some authors 
have even published several papers on heavy-light systems using the same Sum Rule but with 
different results. The difference enters the sum rule through the Wilson coefficients. During this 
work special attention was payed to the mixed condensate coefficient. A short review about the 
different results will now be given. In the first calculation of the J = channel the result was 
incorrect [4] and was corrected in the publication [12]. The sum Rule in the publication [12] is 
today concerned as the correct one in the J = channel. Unfortunately Narison published an 
additional version of the mixed condensate coefficient in the J = which is incorrect [13]. He 
subsequently improves his results step by step until he arrives at the correct one in [18]. There 
have been additional publications in which the Wilson coefficients do not agree with other authors. 
In most cases the quark and gluon condensate coefficients are correct, but for higher dimensions 
the results differ. In the J = 1 channel the situation is similar. This problem was also recognized 
by Ki^ipfer et al. (see [41]). Narison's result in [15] differs from Hayashigaki's [20] in more than 
the mixed condensate terms. The correct result in the J = and J = 1 channel are given in [38] 
and [39]. However, most authors who recognize the ambiguities state that the difference does not 
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change the results significantly [18]. The OPE used for further calculations is given by 
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for the vector current (see [38, 39]). The OPE for the axial partners can be derived from the 
OPEs above by exchanging m to -m. Thus the OPE of the scalar particels can be derived from 
the one of the pseudo-scalar particels and the same holds for the axial-vector and vector parti- 
cles. The terms that are changed by this procedure are the ones which contain odd powers in the 
mass. An inspection of at the OPEs above exhibits that these are terms which contain quarks in 
the condensates. Thus chiral symmetry breaking manifests itself through the quark condensates. 
For vanishing quark condensates the OPEs for the axial partners would be identical. Hence the 
interdoublet splitting is due to these terms. Another important aspect is the mass due to breaking 
of the heavy quark symmetry. Here the intradoublet mass splitting is mainly due to the mixed 
condensate term. In comparison with the finite mass corrections of section 10.1.4 this is quite ob- 
vious because the operator of the mixed condensate is nearly the chromomagnetic interaction term 
which is responsible the biggest contribution to the finite mass corrections and these corrections 
give the mass splittings between the J = and J = 1 states. 

The Wilson coefficient of the unit operator 1 is written down as a dispersion integral. Though 
this notation seems to be quite long winded it is very handy. In many cases the imaginary part 
of the perturbative expression of the current correlator is much simpler than the full expression of 
the current correlator. Moreover if the OPE is applied to a QSR calculation and the continuum 
contribution has to be taken into account it turns out that the notation quoted above implements 
the continuum effects very naturally. The reason is that on the phenomenological side of the OPE 
the continuum of the spectral function is approximated by the imaginary part of the perturbative 
piece of the 2-point correlator, which is included in a dispersion relation. Hence, the imaginary part 
of the perturbative piece of the two point correlator appears both on the left and right hand side of 
a sum rule and is on both sides embedded in a dispersion integral. Thus it can be eliminated from 
the phenomenological side by a simple subtraction and be incorporated by the perturbative Wilson 
coefficient on the theoretical side. In conclusion only the limits of integration in the perturbative 
Wilson coefficient change as the continuum contribution is accounted for. 

The perturbative part of the current correlator for the currents given in (263), with the mass of light 
quarks set to zero, is given in [35]. In that paper a Mathematica package is presented that contains 
the perturbative part of the correlators. The use of the second order corrections is necessary for 
the scalar and pseudo-scalar channel if large q 2 are considered because the imaginary part of the 
correlator becomes negative for large q 2 if only the first order corrections are used. This problem 
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Figure 20: The correlator for the scalar current of a heavy-light system in the approximation 
mught = to first (dashed line) and second order in as (full line). The expressions are nearly 
identical for q 2 < 1000 GeV 2 . Above q 2 = 1000 GeV 2 the curves differ strongly from each other 
until the point is reached where the first order approximation becomes negative. In the case of a 
vector current, the correlators in first and second order remain nearly equal to each other over the 
whole q 2 range. 



is remedied when the second order corrections are included (see figure 20). The requirement of 
positivity for the imaginary part of the correlator stems from the connection between the spectral 
function and the imaginary part of the correlator. The spectral function has to be positive and is 
up to a constant and a power of q 2 equal to the imaginary part of the correlator. Hence, a sign 
switch is forbidden. In fact for small q 2 the expression for the correlator in first and second order 
of ag are nearly equal, while they differ strongly with growing q 2 . This is a beautiful example for 
a calculation where higher order corrections in the coupling constant are necessary in order to get 
a valid approximation. 

In the OPEs (264) and (265) the quark mass effects, as outlined in section 4.7, are taken into ac- 
count. It is a remarkable coincidence that in the limit of vanishing light quark masses the coefficient 
of the gluon and the quark condensate nearly coincide. 

10.2.1 Analysis with moment sum rules 

As the starting point for the analysis of such sum rules, moment sum rules and the narrow resonance 
plus continuum ansatz are used to evaluate the QSR for heavy-light systems. The qualitative result 
of such an analysis concerns the convergence of the sum rule. For D-mesons in the J p = 0~ , 1~ the 
sum rule converges, although the stability of the curves with respect to a shift in the parameters 
is weak. On the other hand in the J p = + , 1 + channels no convergence is seen (see figure 21). 
Hence, there is something wrong, with the Sum Rules in the positive parity channels. The moment 
sum rules in the channel with negative parity is also excluded, but why ? In the example at hand 
the OPEs (264) and (265) with the condensate values of section 5 have been used. The mass of the 
c-quark and the threshold value are the only free parameters. The result should be independent 
from q 2 . The fit of the free parameters should yield nearly equal values independent of the isospin 
or strangeness of the system. Unfortunately this is not the case, the difference in the parameters 
are too large. The plots in figure 22 show the masses for the negative parity states which are 
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Figure 21: Results of the moment sum rules for the J p = + D-mesons. The full line with 
diamonds corresponds to the set q 2 = —lGeV 2 ,m = 1.3GeV,t c = 5QGeV 2 while the other curve 
corresponds to the set q 2 = —5GeV 2 ,m = 1.3GeV, t c — 50GeV 2 . The condensate values are those 
from section 5. No convergence is found for theses sum rules. 



fitted to reproduce the right result. The criteria for the right result was to hit a value which lies 
around the right value for D-mesons with isospin and D s — mesons with strangeness. This should 
lead to nearly equal c-quark masses, a threshold value t c which is close for the two cases and a 
independence from q 2 in a small errorbar window. The threshold requirement is satisfied, but 
all others are not. In addition the plateau for the q 2 = — 1 sum rule is to small to be reliable. 
This does not exclude the sum rule but it diminishes the trust in the reliability of the sum rule 
drastically, without further improvements the sum rule ca not be used. In the J = 1 channels a 
similar behavior is found. 

In the case of B-mesons an analogous behavior is seen, but the convergence in the J p = CP, 1~ 
channel is much better. This is manifested in a wider plateau than in the D-meson case and in a 
better correlation of the parameters. Therefore the moment sum rules for the B-mesons can and 
have been used. At this stage a warning is necessary. In his paper [14] from 1988 Narison calculated 
the mass splittings of the B-meson hypermultiplett with J p = + , CP, 1 + , 1~ of B-mesons with the 
sum rules (264) and (265) and concluded from that the existence of B-mesons with J p = + , 1 + . 
From todays view point that conclusion has to be contradicted. The masses of the positive parity 
states have not been calculated directly, but by Narison's double moment sum rules which calculate 
the quotient of the masses which belong to the states J p — + ,CP or J p = 1 + , 1~. Hence the 
mas splittings of the states J p = + ,CP and the splitting of the states J p = 1 + , 1~ can be 
calculated by his method, but only the masses of the states J p = CP, 1~ can be calculated directly 
by moment sum rules. However no moment sum rule calculation of the D-meson masses has ever 
been published. The only published moment sum rule calculation concerning D-mesons was found 
in the book [15], where an estimate for the decay constants of D-mesons based on moment sum 
rules is given. On the other hand, the calculation of the B-meson masses and decay constants by 
moment sum rules was published separately in a paper (see [15] and references herein). The reason 
for this is probably the bad convergence of the D-meson moment sum rules which have been used. 
There exist many methods to fix the parameters which occur in the moment sum rules. The one 
chosen in this thesis is based on measurements of masses. The decay constants are no directly 
observable quantities, but the masses are. Hence, it can be tested if the sum rule reproduces the 
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Figure 22: Results from the moment sum rules for the J p = CP D-meson with isospin. The full line 
with diamonds corresponds to the set q 2 = — lGeV 2 ,m = 1.3GeV,t c — BOGeV 2 while the other 
curve corresponds to the set q 2 = —5GeV 2 ,m — 1.3GeV,t c = 50GeV 2 . The condensate values are 
those from section 5. The horizontal line gives the experimental mass of CP heavy-light systems 
with isospin. 



experimental mass spectroscopy. If it reproduces the mass spectroscopy and fullfills the stability 
criteria it can be used for the calculation of other quantities like the decay constants. In the case 
of the D-meson moment sum rule it did not prove to be reliable. Thus it can not be used for the 
calculation of the decay constants. 

The moment sum rule in the J p = CP channel is given by: 
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Calculations of D-meson masses based on the OPE given in (264) and (265) using moment sum 
rules are plagued by too many shortcomings, which prohibit high precision analyzes. Thus only 
rough estimates can be obtained by moment sum rules. The problem can be remedied by using 
Borcl sum rules. The discussion of the results will follow in the next section. 



10.2.2 Analysis with Borel sum rules 

The first step to obtain the Borel transformed sum rule is to transform the phenomenological part of 
the sum rule. This leads to a change of the integral kernel from a fraction to exponential function as 
shown in section 6. The imaginary part of the correlator stays untouched by transformation. This 
is the simplest part of the transformation. The transformation of the OPE is in many cases much 
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more work. Here every Wilson coefficient has to be transformed. If the perturbative coefficient 
is given in the form of a dispersion integral at least this one is simple to transform. All other 
coefficients need explicit treatment. The results for the OPEs (264) and (265) are shown below. 
Here are the Borel transformed Wilson coefficients for the pseudoscalar correlator. 



1. Borel transformation of Cq q : 
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Here are the Borel transformed Wilson coefficients for the vector correlator. 
1. Borel transformation of Cq q : 
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2. Borel transformation of Cmixed : 
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This leads to the OPE 
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for the pseudo scalar correlator and 
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for the OPE of the vector correlator. After the transformation the situation for the positive parity 
states changed drastically. The sum rules now converge. Moreover the parameters of all four sum 
rules are now correlated in an acceptable way and the continuum threshold is much smaller. This 
is a welcome effect because the threshold should lie close to the first radial excitation which was 
not the case in the moment sum rule computations. The original analysis of the J = states in 
the strange D-meson channel was done by Narison (see [19] and references herein). Hayashigaki 
performed an analysis of all four states [20]. 

The results of these papers will now be discussed in order to serve as a basis for a further analysis. 
Narison's analysis can be classified as a conservative one. He uses the relative freedom in the 
choice of the continuum threshold t c and the c-quark pole mass m c to obtain agreement with the 
measurements in the D s channel for the J — states. The parameters of his choice are given by 
t c = (7.5 ± 1.5)GeV 2 = (2.725 ± 0.275)Ge^ and m c = lAQGeV. Although the measurements and 
the analysis agreed, the analysis had a huge problem. The Borel curves for the D-meson states 
are expected to have a plateau, as classical mesons like the p had [1], but this was not the case. 
Narison's curves have all been hyperbola shape like. 

Hayashigaki performed an unconservative analysis. He uses a larger value for the charm quark mass 
than Narison and uses another interval for the continuum threshold t c . His analysis reproduces 
the masses of all particles from the hypermultiplet except the + states. For the + states he 
obtained masses larger than the experimental ones. This holds both for the channel with isospin 
and strangeness. He concludes that the £> s (2317) is therefore a four quark state, while he claims 
no clear conclusion for the + state in the channel with isospin. Narison made a comparison 
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of his and Hayashigaki's work in his paper [18], and argues against Hayashigaki. The difference 
between the two publications is probably due to a difference in the perturbative charm quark pole 
mass. Hayashigaki use a perturbative charm quark pole mass of lA6GeV, while Narison used 
1.3GeV. Another problem recognized during the work on this thesis considers the sign of the 
mixed condensate coefficient in the J — 1 channel. In Narison's book [15] he claims a negative sign 
in the case of a vector current, while Hayashigaki has a positive one [20]. Hence, Hayashigaki's 
calculations in the vector channel are in question, while Narison is out of the game because he 
never performed calculations for D-mesons in the vector channel. He did it for B-mesons, which 
are now also in question. One thing holds for Hayashigaki's calculations in all cases, except of the 
1 + case, the curves are hyperbolas and no plateau is visible. Narison just calculated the J = 
channel and there the same phenomenon occurred. The reason for this phenomenon is unclear. 
There are several possibilities 

1. The OPE could be to inexact. The terms could yield too crude approximations or there may 
be terms missing. 

2. The ansatz for the spectral function could be wrong. 

However, there must be something in the sum rules that is right. The likelyhood to write down 
a sum rule which by chance produces results that agree with the measurements is to small to 
seriously consider that possibility. Anyhow it is reasonable to search for possibilities to improve 
the convergence of the sum rule. 

10.2.3 Testing new possibilities to improve the Borel sum rules for the D-meson 
hypermultiplet 

Recent measurements changed the picture of the spectral function owing to heavy-light systems. 
The narrow resonance ansatz is given by a single (^-resonance and a continuum which is a step 
function. The distance between the resonance and the threshold of the continuum lies in the re- 
gion of 0.5 to 1 GcV in the calculations performed by Narison and Hayashigaki. This distribution 
is justified from the J/ip and p spectral function where the spectral functions are measured and 
the assumption holds. During this work no reason for this approximation in other channels was 
found. However, in the axial vector channel of D-mesons there are two cases where the narrow 
resonance plus continuum approximation does not hold. In those spectral functions two states lie 
very close to each other the ground state and a second one, which may be a radial excitation. The 
channels are the s- and u-quark channel of 1+ D-mesons. 

In the s-quark case the resonances are of nearly equal width and mass, while the u-quark case 
consits of a broad and a narrow state (see table 6). The mass splitting is below O.lGeV. Therefore 





1 + 


1 + 


1 + 


V-i) 


m (in GcV) 


2.425 


2.427 


2.459 


2.535 


T (in McV) 


58 


384 


5.5 


2.3 



Table 6: Data of the spectra as they are given by the particle data group (2005). The quantum 
numbers of the heavy particles in the spectral function are not yet sure, the errors are omitted. 

the narrow resonance plus continuum ansatz does not seem to be reasonable. The continuum has 
to start close to the second resonance and if the calculations are performed with Aliev type sum 
rules no convergence is seen with such a small threshold value. Hence, it is reasonable to introduce 
a second resonance and to work with two resonances and a continuum. Moreover, the second 
resonance in the u-quark channel is not narrow but very broad. Therefore the narrow resonance 
approximation is not reasonable any more. As an alternative a Breit-Wigner function is suitable. 
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J (I,S) 


°u 


°(*,o) 


m (in GcV) 


2.255 


2.389 


T (in MeV) 


360 


10 



Table 7: Data of the spectra as they were given in [40] for the + channel with isospin 



Forced by such measurements the question arises if there may be a similar spectrum in the re- 
maining channels + ,0 _ and 1~. Unfortunately no measurements were found, but a theoretical 
prediction. Recently a new method for the calculation of exited states of heavy-light systems has 
been developed by M.F.M.Lutz. Based on this method J.Hofmann calculated the spectral function 
for the J p = + , 1 + states with isospin or strangeness. The result exhibit something in common 
with the mcasurcred spectral functions. Hofmann calculated for six spectral functions the part of 
the spectral function around the ground state. He did it for three channels, one with isospin, one 
with strangeness and one with anti-strangeness. In the channel with isospin he published results 
which can be compared with the measurements. In his publication [40] the axial-vector channel 
with isospin is found which corresponds to the u-quark channel. The calculations which where 
published before the measurements where done agreed qualitatively with the data. Hence, the 
scheme seems to be valid in the positive parity doublet of the D-mesons with isospin. Moreover, 
he calculated further spectral functions. In the + channel with isospin he predicts a spectral 
function which again has two resonance close to each other where one is broad and the other is 
narrow. Therefore a check of the spectral functions with QSR is reasonable for several reasons. 
The quantum number assignment for the measurered spectralfunction can be checked. Moreover, 
the model with which Lutz and Hofmann work does not have quarks and gluons but hadrons as 
basic degrees of freedom. QSR offer a method with which can be checked what QCD says to the 
predictions which have been made. Hence, a test of the method can be performed with such a 
calculation. Therefore the + channel with isospin is checked (see the data in table 7). In order to 
perform the calculations the usual method has to be modified, the parameterization of the spectral 
function must be modified. In the case of the s-quark channel two narrow resonances can be used 
and in the u-quark channel a narrow resonance and a Breit-Wigner resonance can be used 

S(s, m, T) = ^ S 2 T (275) 

7r (s — m z y + s(I y 

where T is the width of the resonance. The continuum is still approximated by the imaginary part 
of the correlator. Thus, the Breit-Wigner curves for the resonances are cut at the threshold t c and 
from there on the imaginary part of the correlator is used. The spectral function for the scalar 
D-meson is than given by 

n (s) = S(s, 2.255GeU, 0.36GeV) + S(s, 2.389GeU, OMGeV) + M + perf (s) (276) 

and plotted in figure23. Due to the Breit-Wigner function the QSR can not be solved for the mass 
corresponding to the broad state. Fortunately, with tools like Mathematica this problem is easily 
manageable. The width of the states will be treated as given by measurements or predictions. Thus 
the narrow resonances can also be approximated by a Breit-Wigner curve without much additional 
work. Figure 23 and equation (276) are not fully correct. The resonances have to carry factors 
which are given by squared matrix elements (0|j|n) 2 . Those matrix have to be inserted into the 
spectral function 

U(s) = (0|j|2.255) 2 S(s, 2.255GeV, OMGeV) + (0|j|2.389) 2 S(s, 2.389Gey, OMGeV) 

+ImU 0+tPert (s). (277) 
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Figure 23: Sketch of the imaginary part of the scalar correlator as predicted by [40] (see the data 
in table 7) . The spectral function contains two resonances a broad one and a narrow one which lie 
close to each other. The continuum is approximated by the imaginary part of the correlator. In 
the case of the axialvector D-mesons the sketch is similar. 



The coefficients can be fitted to reproduce the mass spectroscopy more exactly the ratio of the 
coefficient. This is due to the way the Borel sum rules are evaluated. As always the spectral function 
enters the Borel transformed dispersion integral. In order to eliminate as much as possible from 
the matrix elements in front of the resonances the logarithmic derivative of the dispersion integral 
is taken 



Hence a ratio between the original expressions and the derivative of this expressions is calculated. 
Thus one of the matrix elements can be factored out leaving a ratio between the matrix elements. 
This reduces the number of a priory unknown quantities from three to two. However, there are 
three new parameters in the QSR two widths, one for each particle and the ratio of the matrix 
elements. In QSR calculations with narrow resonances the widths have of course been absent and 
if only one resonance is considered only one matrix element enters the spectral function and this 
is cancelled from the evaluation of the sum rule due to the ratio following from the logarithmic 
derivative. 

Thus the phcnomcnological part of the sum rules is determined. As the theoretical part the OPEs 
of (264) and (265) are used. Therefore a two quark dominance of the states is assumed. The 
expression for the dispersion integral and the OPE are equal. Hence, the logarithmic derivative of 
the dispersion integral and the OPE are equal. This equation can be used to determine properties 
of the spectral function, assuming that the condensates are all known sufficiently exact. Overall 
six free parameters enter the sum rule. 

1. The mass of the first resonance M\. 

2. The mass of the second resonance M 2 . 

3. The width of the first resonance T\. 

4. The width of the second resonance T 2 . 




(278) 
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5. The ratio of the matrix elements r — |o|j|^255) • 

6. The values of the continuum threshold t c . 

Generally it is possible to measure everything except of t c and the matrix elements r. These 
quantities can be measured only in special cases. t c can be measured only for vector currents. A 
particle for which the matrix elements can be measured is the tt meson, where it can be determined 
from decays to a high accuracy. The goal is to check if its possible to fit r and t c so that the sum 
rule reproduces the spectral function. Depending on whether this is possible or not statements 
concerning the physics can be made. 

Before the results are presented a short review of the evaluation of Borel sum rules is given. The 
dispersion relation where on the left side the phenomenological and on the right side the theoret- 
ical part of the sum rules are found is evaluated in order to give the plot of the meson mass in 
dependency of the Borel parameter M. There are two borders for the Borel parameter. The first 
is given by the theoretical part, for small M the OPE is not valid. The second is given by the 
phenomenological part, for large M the continuum contribution gets dominant, but then the reso- 
nance properties can no longer be determined from the sum rule. The Borel window is determined 
to be a domain where the OPE is valid and the resonance dominates the phenomenological part 
of the sum rules. Unfortunately, it is not sure that a Borel window exists. In the case of small M 
the Borel window is fixed, by the requirement that the d = 6 terms contribute less than 10 — 20% 
to the OPE for big M the Borel window is fixed by the requirement that the contribution of the 
resonance to the dispersion integral is bigger than the contribution of the continuum 

/ dslmll(s)e~^ < / ds Imli{s)e~^ 1 . (279) 
Jo Jt c 

Hence the upper border of the Borel window depends on t c , while the lower does not. The Borel 
window determines the regime where the sum rules are reliable in order to determine the meson 
mass. The larger the Borel Window is the more reliable is the corresponding determination of 
resonance properties. 

As it was already discussed above the parameters of the sum rules have to be fixed. Therefore, 
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2.255 ±0.01 
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2.535 ±0.1 


T (in McV) 
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t c (in GeV) 
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t c (in GeV) 


7.5 ±0.3 


7.5 ±0.3 






r (in GeV) 


0.3 ±0.25 


0.3 ±0.25 






m c (in GeV) 


1.2 ±0.01 


1.2 ±0.01 







Table 8: Results from the Sum Rule analysis. Everything except of the widths has been determined 
from the Sum Rules. 

experimental data on the particles are used. The parameters of the Sum Rule are adjusted to re- 
produce the experimental data. Moreover, the Sum Rule has to fulfill additional criteria in order to 
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be reliable. The data which are used are the masses of the particles. The parameters are adjusted 
in order to shift the plateau of the corresponding Borcl curve to the mass of the corresponding 
particle. In addition, the Borel curve is required to be stable against shifts in the parameters. This 
corresponds to an interval in which the parameters can be varied without changing the position of 
the plateau significantly. Finally the plateau has to be located in the Borel window. A large Borel 
window is of course much more reliable than a small Borel window. As the mass of the charm 
quark the running mass was used. The width of the particles was held fixed. 

The intervals for the parameters are given in table 8. In these intervals it is possible to reproduce 
the masses of the corresponding particles. Hence, they can be used to determine the corresponding 
decay constants. After that the Borel window can be determined. Unfortunately, the Borel curve 
in the channel with strangeness shows no plateau in the plots for the decay constants. Due to this 
only one interpretation is possible, the spectral function can not be described by the QSR which 
was used. The other channels have a plateau in the plot for the decay constant, although they do 
not have a plateau in the plots for the masses. The plots for the masses always have the shape of a 
hyperbola. The decay constants agree with earlier calculations. The minimal range for the Borcl 
window for every plot is t = 0.3 — 1.0. Thus, all plateaus lie inside the Borel window. However, also 
the Sum Rules for these states show inconsistencies. The intervals for the parameters are small, 
far from the ideal case. The curves for the masses look like a parabola and not like a plateau, 
some representive curves are shown in figure 24 to 29. Due to these inconsistencies the sum rules 
are also regarded as unreliable. To be more exact, the Sum Rule in the channel with strangeness 
is less reliable than the other Sum Rules. 

The results raise many questions. Even the improved spectral function did not improve the relia- 
bility of the Sum Rules. Moreover, their reliability even got worser. Without the improved spectral 
function the Sum Rule in the channel with strangeness is more reliable than it is here. If the values 
of the decay constants are correct or not can only be proven by experiments. Unfortunately, no 
data from which the decay constants can be extracted are available. The source which is supposed 
to produce the correct result is lattice QCD. The results in these calculations are slightly lower 
than the lattice results, but in agreement with earlier Sum Rule calculations. At this point further 
discussions are not reasonable only further investigations could help to solve the problems. The 
conclusions will offer a possibility to find such a solution. 
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Figure 24: Example for the Borelcurves in the J p = + with isospin. On the left the state with 
m=2.255 GeV and on the right the state with m=2.389 GeV is found. The lines correspond to the 
following parameters: t c = 8.25 for all lines, full line r=4, dashed line r=5 and dotted dashed line 
r=6. 
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Figure 25: Example for the Borelcurves in the J p = + with isospin for the state with m=2.255 
GeV. The lines correspond to the following parameters: t c = 8.25 for all lines, full line r=4, dashed 
line r=5 and dotted dashed line r=6. 
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Figure 26: Example for the Borelcurves in the J p = 1 + with isospin. On the left the state with 
m=2.425 GeV and on the right the state with m=2.427 GeV is found. The lines correspond to the 
following parameters: t c — 7.5 for all lines, full line r=0.1, dashed line r=0.2 and dotted dashed 
line r=0.3. 
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Figure 27: Example for the Borelcurves in the J p = 1 + with isospin for the state with m=2.425 
GeV. The lines correspond to the following parameters: t c = 7.5 for all lines, full line r=0.1, dashed 
line r=0.2 and dotted dashed line r=0.3. 
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Figure 28: Example for the Borelcurves in the J p = 1 + with strangeness. On the left the state 
with m=2.459 GeV and on the right the state with m=2.535 GeV is found. The lines correspond 
to the following parameters: t c — 9.75 for all lines, full line r=l, dashed line r=2 and dotted dashed 
line r=3. 
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Figure 29: Example for the Borelcurves in the J p = 1 + with strangeness for the state with 
m=2.459 GeV. The lines correspond to the following parameters: t c = 9.75 for all lines, full line 
r=l, dashed line r=2 and dotted dashed line r=3. 
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11 Conclusions 



The calculations made in this thesis use an improved ansatz for the spectral function. In most 
of the Sum Rule applications the spectral function is approximated by a single resonance and 
a continuum. In this work a spectral function with two resonances and a continuum is used. 
According to the QCD Sum Rule frame work the calculations are improved when a more realistic 
spectral function is used. Due to the second resonance an additional parameter enters the theory. 
The properties of the resonances are known from experiment or other theoretical approaches. 
Hence, they are used to fix the parameters of the theory. On that basis the decay constants of the 
corresponding hadrons are calculated (see table 8). The results of these calculations and the shape 
of the corresponding Borel curves are used to extract statements concerning the analyzed hadrons. 
The hadrons which are addressed are D-mesons. D-mesons are believed to have two valence quarks, 
a charm quark together with a up, down, or strange quark. Therefore, the Sum Rules used here are 
based on a two quark structure. In particular three systems arc analyzed. A system with J p = + 
with charm and isospin two systems with J p = 1 + one with charm and isospin and the other with 
charm and strangeness. The results are of quantitative and qualitative nature and can be divided 
into two groups. 

In the first group the D-mesons with J p = + and isospin and the D-meson with J p = 1+ and 
strangeness are contained. Although, the Sum Rules for these particles can reproduce the data, 
they are not satisfactory. The OPE does not seem to reproduce the spectral function reliably. 
Although it can not be excluded that an OPE which contains more terms would do it. The 
quantitative results for the D-mesons in the first group are the decay constants in the two quark 
picture. The results for the decay constants of the resonances have again to be split into two 
groups. The first one is the group of the resonances which are the lightest ones. The ground 
state so to say. Their values agree within error bars with earlier calculations. The second group 
is the group of resonances which are higher in mass. No earlier calculations or estimates on their 
value has been found during this work. These properties have most probably never been calculated 
before. For these states the decay constants are always higher than for the ground state particles. 
In the second group the D-mesons with J p = 1 + and isospin is contained. In that channel the OPE 
can not reproduce the spectral function. In consideration of the experience gained during the work 
on this thesis even an OPE with more terms should not change this situation. The calculation of 
decay constants was impossible. The sum rule in that channel did not show any sign of a plateau 
for the decay constant. This point was discussed in section 7. 

However, not long ago some papers on D-mesons and QCD Sum Rules appeared which implement 
a four quark structure of D-mesons [33, 34]. In these publications many problems with which 
the two quark versions are plagued with are absent. Hence, from the viewpoint of the current 
investigation four quark structures have to be taken into account in QSR calculations in the hope 
of maintaining Sum Rules which behave more reasonable. The minimal recommendation which 
can be extracted from this thesis is to consider Sum Rules for D-meson systems which implement 
a two and a four quark structure of D-Mesons in order to extract the mixing angle between those 
structures. This means that it is necessary to consider four-quark structures in addition to the 
standard two-quark ones. 
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